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ABSTRACT 

The SDSS-III/APOGEE survey operated from 2011-2014 using the APOGEE spectrograph, which 
collects high-resolution (R ~ 22, 500), near-IR (1.51 — 1.70/xm) spectra with a multiplexing (300 fiber- 
fed objects) capability. We describe the survey data products that are publicly available, which include 
catalogs with radial velocity, stellar parameters, and 15 elemental abundances for over 150,000 stars, 
as well as the more than 500,000 spectra from which these quantities are derived. Calibration relations 
for the stellar parameters (T e ff, log#, [M/H], [a/M]) and abundances (C, N, O, Na, Mg, Al, Si, S, K, 

Ca, Ti, V, Mn, Fe, Ni) are presented and discussed. The internal scatter of the abundances within 
clusters indicates that abundance precision is generally between 0.05 and 0.09 dex across a broad 
temperature range; within more limited ranges and at high S/N, it is smaller for some elemental 
abundances. We assess the accuracy of the abundances using comparison of mean cluster metallicities 
with literature values, APOGEE observations of the solar spectrum and of Arcturus, comparison 
of individual star abundances with other measurements, and consideration of the locus of derived 
parameters and abundances of the entire sample, and find that it is challenging to determine the 
absolute abundance scale; external accuracy may be good to 0.1-0.2 dex. Uncertainties may be larger 
at cooler temperatures (T e ff < 4000K). Access to the public data release and data products is described, 
and some guidance for using the data products is provided. 
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1. INTRODUCTION 

The third phase of the Sloan Digital Sky Sur¬ 
vey (SDSS-III) included the Apache Point Observatory 
Galactic Evolution Experiment (APOGEE) as one of its 
key components. APOGEE made use of a new, multi¬ 
object, near-IR spectrograph to obtain spectra of more 
than 150,000 stars across portions of the Milky Way vis¬ 
ible from the Apache Point Observatory (APO). The 
main goal of the APOGEE survey is to obtain high- 
resolution spectra to map out the kinematical and chem¬ 
ical structure of the Milky Way. 

The overal l design and goals of the SDSS-III survey are 
described i n Eiscnstein et al. (2011) and of APOGEE in 
iVlajewski et al. J 2015). The survey had a period of com- 
missionmg in spring 2011 and operated in survey-mode 
from fall 2011 until summer 2014. Data from the first 
year of opera tion was released as part of the SDSS DR10 
(Ahn et al. 2014) and included measurements of basic 
stellar parameters as well as overall metal and a-element 
abundances. An SDSS-III internal data release (DR11) 
included the first two years of survey data with the same 
analysis as wa s used for DRIP. The final SDSS-III data 
release (DR12; Alam et al.|2015) includes the full 3 years 
of data from the SDSS III/APOGEE survey. For DR12, 
all of the data have been reprocessed, and the parame¬ 
ter and abundance analysis have been significantly up¬ 
dated: APOGEE DR12 includes measurements of indi¬ 
vidual chemical abundances for 15 different elements. 

The operation of the APOGEE instrument will con¬ 
tinue in the APOGEE-2 survey, which is part of the 
fourth phase of the Sloan Digital Sky Survey (SDSS-IV) 
that plans to operate from 2014 to 2020. APOGEE-2 
will consist of additional observations in the northern 
hemisphere from APO but will also include observations 
of the southern Milky Way that will be obtained with 
the 2.5m du Pont telescope at Las Campanas Observa¬ 
tory. Future data releases for this survey may include 
reprocessing of APOGEE-1 data. 

This paper provides a general overview of the data 
provided by the SDSSTIIM.POGEE survey. )J2] briefly 
summarizes the survey. |3] and [|4] summarize now the 
reduction and spectroscopic analys is is performed (but 
described in greater deta il in Nidever et al. 2015| and 
Garcia Perez et al. 2015) and introduce the reader to 
the associated data products. j|5] discusses the calibra¬ 
tion and validation of the derived - ! stellar parameters and 
abundances. Options for accessing the data are pre¬ 
sented in f(6] iJT] describes the APOGEE catalog data and 
presents some important considerations for users of these 
data, and (j8]has a similar presentation for the spectra. 

2. THE SDSS-III/APOGEE SURVEY 
2.1. The APOGEE instrument 

The APOGEE survey utilizes a fiber spectrograph that 
records 300 spectra simultaneously. The spectrograph 
uses a volume phase holographic (VPH) grating and cam¬ 
era optics to deliver spectra in the H band between 1.51 
and 1.70 /im at a spectral resolution of R ~ 22, 500. The 
spectra are recorded on three Hawaii-2 RG detectors with 
small gaps in spectral coverage between the detectors. 
The instrument was constructed largely at t he Un i ver- 
sity of Virginia and is described in detail in [Wilson et ah 
((2015 ). 


The APOGEE spectrograph was shipped to Apache 
Point Observatory in the spring of 2011 and was com¬ 
missioned during the period April-June of 2011. During 
the commissioning period, it was recognized that there 
was some astigmatism within the instrument, so modifi¬ 
cations were made during summer 2011 to address this. 
Because of the optical issues and modifications after com¬ 
missioning, the commissioning data have lower spectral 
resolution that does not meet the strict survey require¬ 
ments, especially in the detector covering the longest 
wavelengths. Commissioning data have been released in 
the APOGEE data releases as these data provide useful 
radial velocity information, but abundance analysis is not 
currently provided for them because of the reduced res¬ 
olution and quality. However, a number of the sources 
observed during commissioning were reobserved during 
the main survey. 

After the summer 2011 modifications, the instrument 
was pumped down and cooled in August 2011 and re¬ 
mained unopened, operating continuously without any 
large changes in pressure or temperature, for the dura¬ 
tion of the SDSS-III/APOGEE survey. Over the period 
of the main survey, the instrument met survey require¬ 
ments and was remarkably stable, with essentially no 
changes in optical and detector performance. 

At the shorter end of the wavelength range, the detec¬ 
tor pixels do not quite adequately sample the spectro¬ 
graph point-source function. To ensure adequate spec¬ 
tral sampling, the instrument includes a dither mecha¬ 
nism for moving the entire detector array, and exposures 
are taken in pairs, with the detectors shifted 0.5 pixels 
between the two observations in a pair. 

One of the detectors (at the short wavelength end) 
used in the SDSS-III/APOGEE survey exhibits strong 
persistence, where a significant fraction of accumul ated 
charge is re leased over a long period of time (Nidever 
et al. 2015). This affects roughly a third of the chip 
and is distributed in such a way that about a third of 
the fibers are affected at all wavelengths on that chip 
(1.51/im-1.58/im), while the other two-thirds of the fibers 
are unaffected. The behavior of the persistence is chal¬ 
lenging to characterize, and its relevance depends both 
on the brightness of a spectrum as well as the history of 
previous exposures. The current reduction analysis does 
not attempt to correct for this; we present some infor¬ 
mation on the i mpa ct on derived stellar parameters and 
abundances in §5.7[ The offending detector was swapped 
out in summer 2014 after the completion of the SDSS- 
III/APOGEE survey. 

2.2. Survey operations 

The primary source of the li ght to the spect rograph 
is the SDSS 2.5m telescope (Gunn et al. 2006), which 
uses a plug-plate system to collect light from a circular 
3-degree diameter field. For the vast majority of the 
observations, 230 of the fibers were placed on science 
targets, 35 fibers on hot star targets used for correction 
of telluric absorption (although these also enable some 
science projects), and 35 fibers on “blank” sky regions to 
enable sky subtraction. 

Standard exposure sequences consist of 8 separate 500s 
exposures, with 4 at each of the two dither positions, in 
an ABBAABBA pattern. The 4000s integration time 
yields a S/N of about 100 per half-resolution element (a 
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typical pixel) on a source with H = 11, although the 
sensitivity varies across the spectral region recorded on 
the detectors, primarily as a result of grating sensitivity 
changes. 


2.3. lm telescope fiber feed 

Partway through the survey, a fiber feed was con¬ 
structed from the nearby NMSU lm telescope for 10 
fibers, which are located in a fixed linear pattern in the 
focal plane of that telescope. This configuration enables 
observations of single objects, with accompanying sky 
spectra, and has been used to provide a sample of mostly 
bright stars, usef ul for both calibration and science, as 
discussed in §5.5[ To use the lm observations to validate 
those obtained with the 2.5m, the spectra from the two 
telescopes must match closely. The input beam from the 
2.5m telescope is f/5, while that of the lm telescope is 
f/6. However, the fiber run from the lm telescope adds 
significant fiber to the standard 2.5m configuration and 
includes several additional couplings, so the focal ratio 
degradation is expected to be a bit worse for the lm 
feed, which goes in the direction of making the beam 
more closely match that obtained with the 2.5m. 



pixel 



560 565 570 575 580 585 

pixel 


Figure 1. Comparison of sky emission lines as taken with the 
2.5m and lm telescopes. 


Figure [l] shows the profile of two bright sky lines, one 
in the long wavelength chip (chip 0) and one in the short 
wavelength chip (chip 2), for several lm and 2.5m ob¬ 
servations. This demonstrates that the line profiles as 
obtained with APOGEE from the two telescopes match 
extremely well. Comparison of spectra obtained for the 
same star from both telescopes also shows the same good 
match, although the different S/N of the spectra make 
detailed comparison difficult. Both sky and stellar spec¬ 
tra suggest that the NMSU lm delivers extremely similar 
spectra to those delivered by the 2.5m telescope. 


2.4. Target selection 

Targ et selection is described in detail in Zasowsk i et al.| 
(20131. The “main” sample is selected using a simple 
color cut with (J — I\) 0 > 0.5, where colors are dered- 
dened us ing a combination of near-IR and mid-IR pho¬ 
tometry ( Majewski et al.|2011 ), but there are also other 
target classes, including calibration and ancilla ry science 
program targets, as discussed in more detail in [Zasowski] 
etak] ([2013). 

For most of the stars, data are collected in multiple vis¬ 
its to allow detection of radial velocity variables. Stan¬ 
dard survey targets are observed in multiple visits until 
a S/N of 100 per half-resolution element is achieved for 
a target of H=12.2. A minimum of three visits is re¬ 
quired for most targets, with restrictions on the cadence 
so that a broad range of RV variables can be detected. 
However, survey software tracks accumulated S/N, and 
additional visits are added as needed to achieve the S/N 
requireme nt. Some fainter targe ts (medium and long co¬ 


horts; see Zasowski et al. 


2013) were observed for more 


than three visits. A subset of survey targets towards the 
Galactic Bulge was observed only in single visits (with a 
brighter limit for sample selection), because of the chal¬ 
lenge associated with observing such targets from Apache 
Point Observatory, where it is always at a high airmass. 

2.5. Scope of APOGEE data releases 


The APOGEE data in DR10 (Ahn et al. 2014) in¬ 
cluded the first year of APOGEE - observations. The 
SDSS/APOGEE DR11 release added the second year of 
data, but was initially only released internally to SDSS 
collaboration members; it is included in the public DR12. 
Calibration of the DR10 and DR11 data is described in 


Meszaros et al. 

(2013). 



The -SDSS DJ 

412 (Aiam et al. 

2015 

) includes all data 


2011 and July 2014. A total of 163,281 stars are in¬ 
cluded. Of these, 151,141 stars were observed during 
the main survey period that provided survey-quality 
spectra. Hot stars for telluric correction account for 
15,343 of the survey-quality targets, and 17,116 of all 
targets (including commissioning data). Of the 135,798 
survey-quality non-telluric sources, 110,087 come from 
the “main-survey” sample, which does not include spe¬ 
cial targets, targets from ancillary science programs, and 
stars observed with the NMSU lm (883 stars). 

3. BASIC DATA REDUCTION AND ASSOCIATED DATA 
PRODUCTS 

The basic d ata reduction procedure is described in 
Nidever et al. (|2015|). We provide a brief synopsis as 
an introduction to the data products that are included 
in the data releases and described here. 

The raw data consist of data cubes for each individ¬ 
ual exposure, where each layer of the cube is a read of 
the three detectors every 10.7 seconds. As a result, the 
cubes are large, with a standard 500s exposure consist¬ 
ing of 47 reads of 3 detectors, with each detector having 
2560x2048 pixels, since a “reference” array (intended to 
help with tracking bias variations) of 512x2048 pixels is 
included. The raw data cu bes are compressed u sing the 
standard FPACK routine (Seaman et al. 2010), but, to 
provide maximum compression, the raw data are con¬ 
verted to a series of differences between adjacent reads 
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before compression. The resulting files are saved as apR 
files. 

The reduction procedure reduces these image cubes 
to individual 2D images. For each read, dark current 
is subtracted using a calibration product constructed 
from multiple dark images. Up-the-ramp sampling is 
then used to derive count rates that are less affected by 
readout noise than a simple doubly-correlated sampling 
would provide. The resulting image is flat fielded using 
a calibration product constructed from multiple flat-field 
images and is stored in an ap2D file. 

Spectra are optimally extracted from the 2D images 
using an empirical PSF derived from an exposure of il¬ 
luminated telescope mirror covers taken after the sci¬ 
ence exposure sequence has been completed. An approx¬ 
imate relative flux correction is also applied. These ex¬ 
tracted, flux-calibrated spectra are saved as aplD files, 
and the 2D models used for the extraction are saved as 
ap2Dmodel files. 

The multiple exposures (aplD files) are then combined 
to make an individual visit spectrum. To do this, the 
relative shifts between the different exposures, arising 
from the motion of the detector dither mechanism used 
to provide better sampling, are determined. Position- 
dependent sky subtraction is performed on individual ex¬ 
posures, as is a correction for telluric absorption. Expo¬ 
sures at different dither positions are paired, matching 


Stellar parameters and individual chemical abundances 
are derived from the combined APOGEE spectra with 
the APOGEE Stellar Parameters and Chemical Abun- 
da nces Pipeline ( ASPCAP), which is described in detail 


Garcia Perez et al. (2015). The basic idea is that a 


grid of synthetic spectra is searched (with interpolation 
in the grid) to find the best match to each observed spec¬ 
trum, and the parameters of the best match are adopted 
as the best parameters for that star. 

4.1. Synthetic spectra 

Central to ASPCAP is the construction of a synthetic 
spectral grid that covers the range of fundamental stellar 
parameters over which main survey stars are expected to 
be found. Given the significance and frequency of molec¬ 
ular features - in particular, from CO, CN, and OH 
in the region of the iJ-band where APOGEE spec¬ 
tra are recorded, the dimensions of this grid include not 
only the usual effective temperature (T e ff), surface grav¬ 
ity (logg), microturbulent velocity ( v m i cro ), and overall 
scaled-solar chemical abundance ([M/H]), but also the 
carbon ([C/M]), nitrogen ([N/M]), and a-element abun¬ 
dances ([a/M], which includes oxygen) relative to the 
solar abundance ratio. The resulting 7D grid comprises 
several million individual synthetic spectra. 

The fidelity of the synthetic spectral library is criti¬ 
cal to the success of the method. The structure of the 


combined into a spectrum with a sampling correspond- 

ATLAS9 code ( 

Kurucz 

1979 

). For DR10/DR11, a grid 

ing to half the original pixel sampling. The combination 

ot atmospheres was usee 

that adopted solar-scaled a bun- 


technique reduces to simple interleaving (after adjusting 
the levels of the two exposures to match) in the case 
of perfect 0.5 pixel dithers but properly handles dithers 
that are imperfect. The combined pairs are then coad¬ 
ded to provide the final visit spectra. An approximate 
absolute flux calibration is made to the spectrum using 
the 2MASS H magnitude. An initial measurement of 
the radial velocity for each star in the visit is made by 
cross correlating each spectrum with the best match in 
a template library, and these are stored against the visit 
spectra in apVisit files. 

All of the visits are then combined after resampling 
them to a common, log-lambda wavelength scale; in all 
data products, the wavelength scale is in vacuum wave¬ 
lengths. Relative radial velocities of each visit are itera¬ 
tively refined by cross-correlation of each visit spectrum 
with the combined spectrum (which should provide a per¬ 
fect template match). The absolute radial velocity scale 
is then set by cross-correlating the combined spectrum 
with the best match in a template library. The result- 


dances (i.e., the carbon, nitrogren, and a abundances 
were varied only in the spectral synthesis, but not in 
the underlying atmospheric model structure). For DR12, 
a self-consistent set of atmospheres were used: carbon 
and a abundances were varied across the grid for the 
atmosphere calculation (the nitrogen abundance has lit¬ 
tle effect on the structure of the atmos phere, so th is 


was n ot varied). Details are given in Meszaros et al. 
(2012|. These c alculations assume p lane-parallel atmo- 
spheres (but see|Zamora et al. 2015 for some discussion 


of the changes resulting from using a spherical geometry) 
in local thermodynamic equilibrium (LTE). 

The synthetic spectra were generated using the LTE- 
1D branch of the ASSeT synthesis code (K oesterke|2009 k 
A custom lin e list was compiled as discussed in Shetrone 
et al. (2015). The atomic and molecular data for the 


line list were taken from various literature sources, but 
the atomic data transition probabilities were adjusted, 
within limits, by matching synthetic spectra for the Sun 


ing combined spectrum, as well as the resampled visit 

lases ( 

Livingston & Wallace 1991 

Hinkle et al.11995 

) as 

spectra, are stored in apStar files. 

descri 

bed in Shetrone et al. 

(2015). 

The solar abundances 


A summary of the spectral data products is presented 
in Table [T1 The DR12 Science Archive Server (http: 
//data. sdss3. org/sas/drl2) makes both raw data, as 
well as various stages of reduced data, available to the 
community, as discussed in fjoj More details on the files 
are presented in m 

The reduction pipeline software is version-controlled in 
the SDSS-III software repository. 


4. SPECTROSCOPIC ANALYSIS OF APOGEE DATA AND 
ASSOCIATED DATA CATALOGS 


assumed for the se adjustments are th ose from Asplundj 


et al. (2009) and Asplund et al. 

DR12, respectively) and those 
eration of a number of sources; the adopted values for 


([2005]) (for DR10/11 and 
c or Arcturus from consid- 


the va r ious A POGEE line lists a re presented inlShetrone 
et al.| (|2015|). As described in |Shetrone et al.| (|2015|), 
the lineiist has evolved, and future development is likely. 
Additional details on the spectral grids and a compari¬ 
son of the ATLAS9/ASSeT synthetic spectra with those 
from other atmospheres and syntheses are discussed in 
Zamora et al. (2015). 

The high-resolution synthetic spectra are then 
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Table 1 

APOGEE Spectral data products 

Name 

Contents 

apR-[abc]-ID8. apz 
ap2D-[abc]-ID8.fits 
aplD-[abc]-ID8.fits 
apCframe-[abc]-ID8 .fits 

ap Visit- APRED. VERS-PLATE-MJD-FIBER.fits 

apStar-A PRED_ VERS-APOGEE_ID.fits 
aspcapStar-RESULTS-VERS-APOGEE JD. fits 

Raw data cubes (2560x2048x47), APOGEE compressed format 

2D (2048x2048) dark subtracted, flat fielded 

ID extracted (2048x300), rough wavelength calibration (not zero-pointed using sky) 

ID spectra with sky subtraction, telluric correction, wavelength zeropoint adjustment 
Individual visit spectra, dither-combined, but otherwise native pixel scale, 
one per object visit 

Spectra resampled to constant A log A, along with combined spectrum, one per object 
Normalized spectra to match those of spectral grids, 
with best-matching grid spectra, one per object 


matched to the APOGEE resolution and pixel sam¬ 
pling. This is accomplished by smoothing the high- 
resolution spectrum with the line-spread function (LSF) 
of APOGEE data. However, the APOGEE LSF varies 
across the detector as a consequence of the optical de¬ 
sign. To some extent, this variation is reduced when 
multiple visits are combined, since stars will generally 
land in a different location on the detector (from a dif¬ 
ferent fiber in the telescope focal plane) in different visits. 
In DR10, a Gaussian LSF with R = 22500 was adopted, 
but we subsequently recognized that the APOGEE LSF 
is not well-represented by a Gaussian. For DR12 we 
have taken a approach of using an LSF made by av¬ 
eraging five different empirically-measured LSFs across 
the detector to make a single library for use with all of 
the SDSS-III/APOGEE spectra. As demonstrated be¬ 
low, this yields results that are validated by observa¬ 
tions of objects of known abundances, but it is an area 
for poten t ial i mpro vem ent in subsequent data releases. 


Garcia Perez et al. (2015) present a more detailed inves- 
tigation of the potehtial impact of assuming an imperfect 
LSF. 

In addition to smoothing the high-resolution synthetic 
spectra with an APOGEE LSF, in DR12 we also include 
broadening to account for macroturbulent velocity in the 
stellar atmosphere. We were motivated to do this not 
only because it is a known physical phenomenon, but 
also because this additional broadening brings our de¬ 
rived surface gravities closer to those suggested by other 
techniques (see below). For DR12, we have adopted a 
simple scheme of broadening with a Gaussian of FWHM 
6 km/s for the m acrotur b ulent veloci ty, motivated by 
typical values from Carney et al. (2008). Adopting a dif¬ 
ferent functional form, along with possible variations in 
macroturbulent velocity as a function of other stellar pa¬ 
rameters, is another area for potential improvement in 
the future. 

The synthetic spectra are calculated relative to a true 
continuum. However, it is challenging, if not impossi¬ 
ble, to identify a true continuum in the observed spectra 
in cool, metal-rich stars. To account for this, both the 
synthetic spectra and the observed spectra are “pseudo- 
continuu m” normalized using the sa me algorithm, as dis¬ 


cussed in Garcia Perez et al. 


(2015). 


The smoothed, resampled, and pseudo-continuum nor¬ 
malized spectra are bundled into a large library. To min¬ 
imize the (substantial) computing resources required for 
such a large grid, we separate the grid into two grids 
overlapping in temperature (which we denote GK and 
F); we find the best match in each grid, and use the 


Table 2 

Synthetic spectral grid parameter values 


gridname Dimension 

Range 

Step size 

GK grid T eff 

3500 to 6500 

250 

log g 

0 to 5 

0.5 

log Vmicro 

-0.301 to 0.903 

0.301 

[M/H] 

-2.5 to 0.5 

0.5 

[C/M] 

-1 to 1 

0.25 

N/M 

-1 to 1 

0.5 


a/ M] 

-1 to 1 

0.25 

F grid T e ff 

5500 to 8000 

250 

log 9 

0 to 5 

0.5 

log Vmicro 

-0.301 to 0.903 

0.301 


M/H] 

-2.5 to 0.5 

0.5 


C/M] 

-1 to 1 

0.25 


N/M] 

-1 to 1 

0.5 


a/M] 

-1 to 1 

0.25 


result that provides the better match. In addition, the 
libraries are compressed using a principal components de¬ 
composition, in which we split the spectra into 30 pieces 
and PCA-compress each of those, keeping only the first 
30 principal components and the amplitudes of these for 
ea ch location in the grid; ad ditional details are provided 


Garcia Perez et al. (2015). 

The range and spacing of the ATLAS9/ASSeT grid 
par ameters us ed for DR12 are given in Ta ble [2] (see also 
Zamora et al. 2015). Garcia Perez et al. (2015) discuss 


why this grid spacing is adequate and how interpolation 
within the grid is performed. 

We note that for SDSS-III/APOGEE, stellar rotation 
is not included in the synthetic spectral grids. For stars 
with rotational velocities (> 5 km/s) that broaden the 
lines significantly beyond the instrument resolution, the 
libraries will not be able to obtain a good match to the 
observed spectra, and the resulting stellar parameters 
can be affected systematically. Stellar rotation is not ex¬ 
pected to be significant for most of the red giants that 
make up the APOGEE sample. However, a significant 
fraction of the dwarfs observed in the survey may have 
detectable rotation, and results for these are suspect. As 
discussed below, we do not provide calibrated parame¬ 
ters for dwarfs, partly because of the lack of treatment 
of rotation. For all survey stars, we provide a flag that in¬ 
dicates whether rotation is likely to be important, based 
on a comparison of the width of the autocorrelation func¬ 
tion of the spectrum relative to the autocorrelation of 
the best-matching RV template. Future improvements 
are planned with the inclusion of rotation when needed. 

We also note that the current analysis is limited to 



































6 


[M/H] 





O 

o 

1 

> 


3 

2 

1 



-2.5 


- 2.0 


-0.5 


2 3 4 5 

logg 


Figure 2. Microturbulence relation, as derived from a 7D run of 
the calibration subsample. 


stars with T e g> 3500 K because the Kurucz grid of 
model atmospheres does not extend to cooler temper¬ 
atures. Subsequent work is planned to extend results to 
cooler temperatures using MARCS model atmosp heres 


and t he Turbospectrum synthesis code (Zamora et al. 


20151. 


4.2. Deriving stellar parameters and abundances 

For each APOGEE spectrum, a best match within 
the synthetic grid is found. To do this, interpolation 
within the grid, as well as an efficient search method, 
are required. This is acco m plishe d using the FOR- 
TRAN95 code FERRE (Allende Prieto et al.||2006plL, 
which yields a best set of stellar parameters. The com¬ 
putational resources to perform the full 7D search for 
> 100,000 stars are substantial (> 50,000 CPU hours). 
Because of this, we use a calibration subsample (see (j5| 
to derive a relation for the microturbulent velocity as a 
function of surface gravity. Figure [2] shows the derived 
microtubulent velocity as a function of surface gravity 
for the calibration subsample, color-coded by metallicity. 
The stars of higher metallicity that make up the bulk of 
the APOGEE sample (but not the calibration sample) 
define a fairly tight sequence. We a fit a linear relation 
to the subset of stars from these data with logg< 3.8 
and [M/H]> —1 and derived the fit that was adopted for 
DR12: 

Vmicro = 2.478 - 0.325 log g, (1) 


which is shown as a solid line in Figure [2] We then use 
this relation to construct a 6D grid, interpolating within 
the 7D grid to the microturbulent velocity appropriate 
for the surface gravity of each grid point. The utility of 
this is that it reduces the search to a 6D one, rather than 
7D, with a corresponding reduction in the time required 
(by a factor of three or more) to analyze the full sample. 

The full DR12 sample is then analyzed using the 6D 
grid to provide the initial set of parameters: T e ff, log g, 
[M/H], [C/M], [N/M], and [a/M]. 

Given the stellar parameters, the spectra are then an¬ 
alyzed for abundances of each of 15 different elements — 


C, N, O, Na, Mg, Al, Si, S, K, Ca, Ti, V, Mn, Fe, and 
Ni — that have features in the APOGEE wavelength 
coverage. For each element, the regions in the spectrum 
that are sensitive to the abundance of that element, and 
less sensitive to the abundance of other elements, are de¬ 
termined from model spectra. The spectral “windows” 
identified in this way are weighted by their relative sen¬ 
sitivities and also by how well a synthetic spectrum with 
the parameters appropriate to Arcuturus match those of 
an observed Arcturus spectrum. Details of this proce¬ 
dure, and the s pecific locations of the der ived windows, 
are discussed in Garcia Perez et al. (2015). 

Given the windows tor each element, each spectrum 
is then run through the FERRE search code separately 
for each element, using the same grid as used for the 
determination of the stellar parameters. However, for 
the individual element abundances, the search is done 
holding the stellar parameters T e g, log g, v m i cro fixed, 
and just varying a single metallicity parameter, so this is 
a one-dimensional search. For C and N, the [C/M] and 
[N/M] dimension is varied, for the a elements (O, Mg, 
Si, S, Ca, and Ti), the [a/M] dimension is varied, while 
for all others (Na, Al, K, V, Mn, Fe, Ni), the [M/H] 
dimension is varied. Apart from C and N, this means 
that more than just the single element abundance being 
sought is being varied, but ASPCAP attempts to account 
for that by using windows with maximal sensitivity to the 
element being derived and minimal sensitivity to other 
abundances. This allows us to avoid having to create 
a series of 15 different 7D grids where pure individual 
element variations are calculated. 

4.3. FERRE output 

The full ASPCAP pipeline produces an initial set of 
six parameters for each star from the initial 6D fit, plus 
15 additional abundance measurements. The raw output 
from the 6D fit are stored in an array called FPARAM, 
while the output from the 15 separate ID “window” fits 
are stored in an array called FELEM. 

FERRE also produces a covariance matrix from the 
initial 6D fit, which we store in a matrix called 
FPARAM_COVAR, and individual abundance uncer¬ 
tainty estimates, which we store in a FELEMJERR ar¬ 
ray. Unfortunately, we have found that the formal errors 
returned by FERRE seem to be significantly underesti¬ 
mated when they are calculated using standard fit tech¬ 
niques (i.e., inverse of the curvature matrix); this is plau¬ 
sibly caused because systematic errors (e.g., in the syn¬ 
thetic calculat ions, LSF-matching, etc.) dominate statis¬ 
tical ones (see Garcia Perez et al.|2015 for more details). 
As a result, we also determine ana store a more empirical 
estimate of uncertainties, as discussed below. 

For each parameter and elemental abundance, we store 
a bitmask flagging possible conditions encountered in the 
process, such as a solution falling near the edge of the 
grid, and consequently deemed to be less reliable. 

The data products are discussed in more detail in £j6] 
The ASPCAP software is version controlled in the SDSiT 
III software repository in the speclib, idlwrap and ferre 
products. 


4.4. Versions and data releases 


1 http: //wire, sdss3.org/svn/repo/apogee/aspcap/ferre, 

http://hebe.as.utexas.edu/ferre 


There are four main stages of APOGEE data 
processing: visit reduction, star combination, ASP- 
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Figure 3. Temperature calibration, derived from low-reddening 
sample with photometric temperatures. 

CAP/FERRE processing, and final calibration. The 
software is set up such that these four stages can be 
run with independent software versions/configurations so 
that, for example, a later stage can be run without hav¬ 
ing to rerun an earlier stage with a new software version. 
The output from the four stages is saved in a series of 
nested subdirectories, as described in S|6] The names of 
the different versions are tied to software version num¬ 
bers of the different products. Table [3] gives the names 
and versions used for the APOGEE data releases. 


5. CALIBRATION AND VALIDATION OF DR12 
PARAMETERS AND ABUNDANCES 

Results for the stellar parameters and abundances have 
been calibrated and validated based on observations of 
a calibration sample. This sample includes stars in the 
Kepler field that have well-characterized effective tem¬ 
peratures and surface gravities, samples of stars in a set 
of globular and open clusters that span a wide range of 
metallicity, and a set of stars with previously-measured 
parameters and abundances. 


5.1. Effective temperatures 

We compare the ASPCAP spectroscopically- 
determined temperatures with those obtained from the 
cali b rated photometric relation of Gonzalez Hernandez] 
& Bonifacio (2009). We chose a sample of giants (logg< 
3.8) in the Kepler field and adopted E(B-V) estimates 


as de termined from the SAGA survey (Casagrande et al. 
2014|. We convert E(B — V ) to E(J — K ) using 


E(J -K) = (0.56 + 0.06 (J - K) 0 ) * E{B - V) 


(Bessell et al. 1998). We use the dereddened J — K color 
to derive the temperature, and adopt the raw ASPCAP 
metallicity for the color-temperature relation. We re¬ 
strict the comparison to stars with ASPCAP metallicity 
[M/H]> — 1. Figure [3] shows the difference between the 
spectroscopic and photometric temperatures for this sub¬ 
sample. This comparison suggests that there is an offset 
in the derived effective temperature in that the raw AS¬ 
PCAP temperature is about 90 K cooler than the pho¬ 
tometric temperatures. 

We derived a calibration based on a linear fit with tem¬ 
perature to these data, which yields: 


Teff.corr = T eff , raw - 0.0034(T eff , raw - 4500) + 87.9 (2) 

although the dependence of the difference on tempera¬ 
ture is negligibly small. 


The calibrated temperature provided for all giants 
in the DR12 release uses the calibration from equa¬ 
tion [2] For giants with effective temperatures out¬ 
side the range of those used to derive the calibration 
(3500 <T e ff< 6000), the correction at the edge of the 
range was adopted, and the CALRANGE_WARN bit has 
been set in the PARAMFLAG for temperature (see <Q. 

An upper limit on the temperature uncertainty is de¬ 
rived from the scatter around this relation, which gives 
a(T e ff) = 91.5 K. This is an upper limit because there 
are uncertainties associated with the photometry. On the 
other hand, they may be some systematic uncertainties 
related with the photometric temperature relation. This 
scatter around the calibration relation was adopted as 
the empirical temperature uncertainty f or stars in DR 12. 

The color- t empe rature relation ofjGonzalez Hernandez| 
& Bonifacio (2009) is also used to nag objects that may 
have poor "stellar parameters. For each object, a temper¬ 
ature is determined using the color-temperature relation, 
the 2MASS J — K color, and the estimate of the redden¬ 
ing used by the targeting pipeline. This is compared with 
the ASPCAP-derived temperature. When these deviate 
by more than 1000 K, the COLORTE_WARN bit is set 
in the ASPCAPFLAG (see Q; if the temperature differ 
by more than 2000 K, the COLORTE_BAD bit is set and 
the ASPCAP results are flagged as bad. Most of these 
are cool dwarfs that are not fit well (in retrospect, we 
might have chosen less conservative values to flag). 

5.2. Surface gravities 

Surface gravity is a key parameter for a number of 
reasons: (1) it is important for estimating distances to 
survey stars, (2) it can be used to identify AGB stars, 
and (3) it can have a significant effect on line strengths. 

The APOGEE survey includes observations of giants 
that have time series photometry with Kepler that allows 
asteroseismic analysis. In principle, such analysis very 
accurate surface gra vities, with an estimated uncertainty 
of << 0.05 in logg ( Hekker et al.||2013 |. 


A calibration relation tor surface gravity has been de¬ 


APOKASC catalog 

Pinsonneault et al. 

2014 

. To cover 

a wide range of prop 
ferrable evolutionary 

erties, we took all si 

ars t 

lat had in- 

states from Stello et al. 

(2013 

I, all 


stars with logg<2, and all of the small number of stars 
in the APOKASC catalog with [M/H] <-l. 

Figure 0] shows the difference between ASPCAP spec¬ 
troscop y Tog a and the asteroseismic log g from P inson-| 
neault et al. ( |2014| ), as a function of surface gravity. It 
is clear that results for stars at higher surface gravity 
(logs > 4) are poor. We expect that this is caused to a 
large extent by the lack of treatment of rotation discussed 
above. 

At lower surface gravities, there are systematic offsets 
between gravities determined for RGB and RC stars; the 
difference between ASPCAP and asteroseismic gravities 
is larger for RC stars, even for RGB/RC pairs at the same 
gravity and metallicity. This has alrea dy been no ted in 
previous APOGEE samples (Pinsonneault et al. 2014 


Meszaros et al. 20131), but the effect is still currently not 
well understood! If one considers just the RGB sam¬ 
ple, it appears that the difference between ASPCAP and 
asteroseismic gravity increases to lower gravities. No de¬ 
pendence on metallicity was found. 















































Table 3 

Software version and product names 


Data release 

APREDJVERS (apogeereduce) 

APSTAR_VERS (apogeereduce) 

ASPCAP.VERS (idlwrap) 

RESULTS.VERS (idlwrap) 

DR10 

r3 (v2_ll) 

s3 (v2_ll) 

a3 (v7) 

v304 (vl5) 

DR11 

r4 (v2_14) 

s4 (v2_15) 

a4 (vl6) 

v402 (-) 

DR12 

r5 (v3_04) 

stars (v3_10) 

125_6d (v22) 

v603 (v27) 


Table 4 

Calibration clusters 


Figure 4. Surface gravity calibration, derived from a sample of 
stars in the Kepler field with gravities from asteroseismic analy¬ 
sis. This analysis also provides information on the evolutionary 
state: red points are hydrogen shell burning (RGB) stars, while 
blue points are helium core burning (RC) stars. 

We have chosen to adopt a surface gravity calibration 
using a linear surface gravity fit to the RGB sample. If 
one assumes that the asteroseismic gravities are correct, 
this would imply that the RC stars in the DR12 sample 
will have calibrated surface gravities that are too large by 
~ 0.2 in log<?. We have chosen this approach rather than 
“splitting the difference” between the RGB and RC sam¬ 
ples because it is possible to distinguish the bulk of the 



relation that we hope will give the best results for RGB 
stars, and recommend treating the RC stars separately. 
The derived calibration relation is given by: 

^g Qcorr = lo§ draw T 0.14 log (/ rau , 0.588. (3) 

For uncertainties in logg, we adopt the observed scat¬ 
ter around this calibration relation, again ignoring the 
RC stars. This gives an empirical uncertainty of 0.11 in 
log g. However, the true error is probably more likely 
to be dominated by systematics, given the Kepler differ¬ 
ences. We have also determined that the surface gravity 
is sensitive to the details of the LSF used in the analysis, 
so this may lead to systematics depending on the exact 
LSF combination of each spectrum. 

5.3. Calibrated parameters and abundance 
determination 

Our procedure to determine elemental abundances uses 
the raw derived stellar parameters when fitting synthetic 
spectra to the windows in the spectra that cover features 
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lull: 


of specific elements. We choose to do this because (1) 
using calibrated parameters would lead to significantly 
poorer overall spectrum fits, and (2) it is possible that the 
raw, “spectroscopic,” parameters account to some extent 
for uncertainties and approximations that are implicit in 
the synthetic spectra, and so using the raw parameters 
provides a more self-consistent set of abundances. 

On the other hand, it is also possible that the raw pa¬ 
rameters are offset for other reasons, for example, if the 
adopted LSF leads to incorrect surface gravities. If this 
were the case, then one might argue that using the cal¬ 
ibrated parameters for individual elemental abundances 
might be more accurate, at least for features where the 
strength depends strongly on surface gravity. Future 
work is planned to investigate this. 

5.4. Chemical Abundances 

Our primary calibration sample for abundances is a set 
of globular and open clusters spanning a wide range of 
metallicity, from M92 and M15 on the metal-poor end to 
NGC 6791 on the metal-rich end. The set of calibration 
clusters is presented in Table [4j 

We have approached calibration in two steps: an in¬ 
ternal calibration where we look for homogeneity of ele¬ 
mental abundances within clusters, and an external cali¬ 
bration where we investigate the agreement of the abun¬ 
dances with previously measured values. 

5.4.1. Internal calibration 
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For the internal calibration, we examine the ASPCAP 
abundances under the assumption that clusters are chem¬ 
ically homogeneous for all elements except carbon and 
nitrogen, which are expected to have variations in giants 
because of potential mixing of CNO-processed material. 
The sample was restricted to more metal-rich clusters be¬ 
cause (1) multiple populations in some globular clusters 
are manifest in elemental abundances (e.g., Al, Mg, Na, 
and O), and our cluster sample is dominated by globu¬ 
lar clusters at [M/H]< -1; (2) at lower metallicity, the 
lines are weaker, so trends are more challenging to dis¬ 
cern; and (3) the bulk of the main APOGEE sample is 
at higher metallicity, so we choose to provide an inter¬ 
nal calibration that is most accurate for the bulk of the 
sample. In practice, we inspected the spread of each in¬ 
dividual element within clusters of various metallicities, 
and chose a lower limit for each element to use for the in¬ 
ternal calibration; this limit was always between [M/H]of 
-1 and -0.6. 

For the internal calibration, we look for agreement 
among the cluster stars, without any constraints on what 
the mean abundances are. Inspection of the raw ASP- 
CAP abundances within each cluster suggests that the 
abundances of some elements appear to have small trends 
with effective temperature. To quantify this effect we fit, 
for each element, for mean abundances of each cluster 
plus a linear temperature dependence of the abundance: 


[Xy/H] = [Xi/H] + S x (Tij - 4500)/1000. (4) 

where i represents the cluster, and j the star within the 
cluster, and Sx the derived slope for element X. Note 
that we used [X/H] for all abundances: for the elements 
for which ASPCAP uses a grid dimension that is relative 
to the overall metal abundance (C, N, and a elements), 
we added the ASPCAP [M/H] to get [X/H]. This made 
it easier to detect potential cluster non-members as those 
stars with abundances significantly different from the 
cluster averages. We attempted to be conservative about 
cluster membership: we starte d with stars selected to be 
probable cluster members (see Zasowski et al. 2013) , but 
then adopted only stars with radial velocities within 5 
km/s of our cluster mean, and then finally rejected stars 
with abundance residuals more than 0.2 dex from the fit. 

Results are shown in Figure [5] where abundances of 
each cluster calibration star are shown relative to the 
derived mean ([X;/H]) of the cluster of which it is a 
member. Points are color-coded by the metallicity of 
their respective cluster. The small temperature trends 
could arise from a number of different factors: the as¬ 
sumption of LTE, issues with atomic data since different 
lines provide information about abundances at different 
temperatures, etc. The data suggest that a linear tem¬ 
perature term is adequate to describe the trends for most 
elements. 

For DR12, we have adopted the derived temperature 
trends for an internal calibration of the abundances for 
all elements except carbon and nitrogen. The slopes of 
these temperature terms (Sx) are given in Table [5] We 
arbitrarily define the abundance of a given element as 
the abundance at T e g= 4500 K. While is it possible that 
the derived abundances of carbon and nitrogen may have 
temperature dependencies, we cannot discern these from 
the cluster data internally. 


Table 5 

Internal temperature 
calibration slopes 


Element slope per 1000 K 


c 


N 


O 

+0.007 

Na 

-0.001 

Mg 

-0.051 

Al 

-0.046 

Si 

-0.121 

S 

+0.060 

K 

-0.046 

Ca 

+0.059 

Ti 

-0.162 

V 

-0.177 

Mn 

-0.206 

Fe 

-0.055 

Ni 

-0.124 

[M/H] 

-0.056 

[a/M] 

-0.028 


We apply the temperature corrections as derived from 
the cluster sample to the DR12 sample. Because the 
clusters only provide stars with 3800 <T e g< 5250K, we 
do not extrapolate the linear relation outside this range. 
Instead, we use the correction at T e ff = 5250K for all 
stars with T e ff > 5250K, and the correction at T e // = 
3800K for all stars with T e ff < 3800K. For those giants 
outside the temperature range of the calibrator, we set 
the CALRANGEJWARN bit in the abundance masks. 

In addition, because the calibration is derived only 
from giants, and because of other issues noted above for 
dwarfs, we do not apply it to stars with logg > 3.8. For 
these higher surface gravity stars, we do not provide cal¬ 
ibrated abundances at all (the values are set to -9999) 
and we set the CALRANGEJ3AD bit in the abundance 
masks. The raw ASPCAP parameters are still provided. 

5.4.2. Empirical abundance precision 

The individual abundances of stars in clusters provide 
an excellent data set to assess the precision of our derived 
abundances, by using the scatter within clusters as an 
estimate of the abundance uncertainties. This scatter is 
a function of both temperature and metallicity, which is 
as expected because lines are weaker at lower metallicity, 
and change strength with temperature. In addition, we 
see that scatter decreases with increasing signal-to-noise. 

To derive empirical uncertainties, we calculated the 
spread of the abundance within each cluster in bins of 
width 250K in temperature and 50 in S/N (putting all 
observations with S/N> 200 into the highest S/N bin). 
To achieve a greater range in S/N, we included abun¬ 
dance measurements for cluster stars derived from in¬ 
dividual visits in addition to those from the combined 
spectra. We combined results for all clusters in bins of 
0.5 dex in metallicity, and then fit a simple functional 
form to the uncertainty in each element: 


ln(crx) = A x + B x (T e f f - 4500)+ 

C x [M/H] + D x ( S/N - 100) (5) 

where we fit for the natural log to ensure that the uncer¬ 
tainty never goes negative (i.e., outside of the range of 
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Figure 5. Internal elemental abundance calibrations. Plots show deviations of individual elements relative to the cluster means as 
a function of temperature; different symbols are for different clusters, as noted in lower right. Lines shows the adopted temperature 
calibrations. Carbon and nitrogen are not shown because we cannot assume homogeneous abundances within clusters. 
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[Fe/H] ASPCAP 

Figure 6. External elemental abundance calibration, 
the calibration cluster data). 

Table [6] gives the coefficients for these empirical er¬ 
ror estimates for all of the abundances. The second-to- 
last column gives the empirical error estimate at T e ff= 
4500K, solar metallicity, and S/N = 100. For all ele¬ 
ments except vanadium, the observed scatter is less than 
0.1 dex, and for many elements, it is below 0.05 dex. 
Generally, the errors increase with increasing tempera¬ 
ture, decreasing metallicity, and decreasing S/N. 

These empirical errors, which are distributed in the 
APOGEE data release (along with the raw FERRE er¬ 
rors), capture the scatter in elemental abundances ob¬ 
served in clusters over a narrow range of effective tem¬ 
perature. Over a broader range of temperature, they do 
not capture the degree to which our simple, metallicity- 
independent, linear temperature calibration relation fits 
the data. Since the observed stars in any given cluster do 
not usually span a large temperature range, we estimate 
the overall uncertainty in abundances, across a broad 
range of temperature, by measuring the scatter around 
the calibration relations of Figure [5j This is presented in 
the final column of Table [bj These “global” uncertain¬ 
ties do not capture the likely dependence of the precision 
on temperature, metallicity, and signal-to-noise. They 
represent a more conservative estimate of the internal 
uncertainties. 

The scatter within clusters indicates that the precision 
of our abundances is less than 0.1 dex for all elements, 
with the possible exception of vanadium. For some abun¬ 
dances, especially if considering the sample over a re¬ 
stricted temperature range and high S/N, the precision 
can be better than 0.05 dex. 

5.4.3. External calibration 

Table El gives the adopted literature values for [Fe/H], 
along with references, for the calibration cluster sample. 
The difference between the derived ASPCAP values and 
the literature values for both the parameter [M/H] and 
the abundance [Fe/H] are shown in Figure [6] We adopt 
the literature [Fe/H] as a reference for our [M/H] values 
because our method allows for variations in a elements, 
carbon, and nitrogen separately; implicit in this compar¬ 
ison is the assumption that the Fe abundance drives the 
overall metallicity in the APOGEE spectra. 

The derived ASPCAP [M/H\ and [Fe/H] are similar 


to the literature values near solar metallicity, but are 
~ 0.2 — 0.3 dex more metal-rich than the literature for 
very metal-poor clusters. This is similar to what wa s 
seen for the DR10 data release (Meszaros et al. 2013). 
We adopted a quadratic relation m metallicity to derive 
calibration relations: 


[M/H] corr = [M/H] raw + 0.0260+ 

0.255[M/H] raw + 0.062 [M/H] 2 raw (6) 

[ Fe/H\ corr = [Fe/H\ raw + 0.0326+ 

0.245 [M/H\ raw - 0M2[M/H] 2 raw (7) 

We have chosen to apply the calibration above to the 
[M/H] parameter, but not to the [Fe/H] abundances, 
given the impossibility of calibrating all of the abun¬ 
dances in a homogeneous fashion. Most users will pre¬ 
fer to work with [X/Fe] by subtracting the [Fe/H] from 
[X/H], and in many cases systematic errors in [Fe/H] 
and [X/H] will partly cancel out. If we were to apply the 
external calibration to [Fe/H], it would have the effect 
of increasing the [X/Fe] ratios of all elements towards 
lower metallicity. Since we do not know the source of the 
[Fe/H] discrepancy at low metallicity, it is difficult to 
judge whether it is likely to affect all elements similarly, 
or whether it is specific to iron. As a result, we choose 
to leave all of the [X/Fe] ratios on the native, consistent, 
ASPCAP scale. 

The main issue addressed by the calibration relation 
is at the metal-poor end. If one were to ignore all clus¬ 
ters with [Fe/H]< —0.5, the remaining metal-rich clus¬ 
ters show little evidence for an offset or a trend with 
metallicity. However, the simple quadratic form of the 
relation adopted to address the issues at [Fe/H]< —0.5 
results in a relation that changes by almost 0.3 dex 
at —0.5 <[Fe/H]< 0.5. As a result, there are non- 
negligible differences between the calibrated and uncali¬ 
brated [M/H] in this regime, and the uncalibrated value 
may be just as good as the calibrated ones for clusters 
with [M/H]> —0.5. 

5.5. The abundance scale 

The ASPCAP abundances are presented as relative to 
hydrogen in the DR12 database in the usual bracket no¬ 
tation: 


[X/H] = log (n x /n H ) - log (n x /n H )o 


( 8 ) 


where n-x and nn are the number density of nuclei of el¬ 
ement X and hydrogen, respectively. But some caveats 
apply. The abundances are not truly differential to the 
values that APOGEE derives for the Sun. Solar abun¬ 


dances are adopted from Asplund et al.] (2005) and used 
for computing model atmospheres (Meszaros et al. |2012 1 
and fine-tuning the line list (Shetrone et al.||2015|). The 


line list used for spectral synthesis includes, in addition 
to laboratory and theoretical transition probabilities and 
damping constants, modifications to provide a compro¬ 
mise match to both the solar spectrum and Arcturus. 
However, only atomic lines were adjusted, and even for 
those, only within limits. 
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Table 6 

Empirical uncertainty relation coefficients and values 


Element 

A 

B 

C 

D a (T eff 

= 4500,[M/H]= 0,S/N= 100) “Global” uncertainty 

c 

-3.350 

0.769 

-0.919 

-0.066 

0.035 


N 

-2.704 

0.291 

-0.591 

-0.078 

0.067 


O 

-3.649 

0.670 

-0.614 

-0.093 

0.026 

0.050 

Na 

-2.352 

-0.002 

-0.915 

-0.263 

0.095 

0.064 

Mg 

-3.537 

0.263 

-0.825 

-0.297 

0.029 

0.053 

Al 

-2.764 

0.471 

-0.868 

-0.162 

0.063 

0.067 

Si 

-3.150 

0.383 

-0.224 

-0.105 

0.043 

0.077 

S 

-3.037 

0.507 

-0.625 

-0.299 

0.048 

0.063 

K 

-2.770 

0.216 

-0.667 

-0.275 

0.063 

0.065 

Ca 

-3.226 

0.284 

-0.879 

-0.429 

0.040 

0.059 

Ti 

-3.186 

0.657 

-0.819 

-0.068 

0.041 

0.072 

V 

-1.608 

0.900 

-0.400 

-0.418 

0.200 

0.088 

Mn 

-3.031 

0.639 

-0.661 

-0.326 

0.048 

0.061 

Fe 

-3.357 

0.098 

-0.303 

-0.071 

0.035 

0.053 

Ni 

-3.153 

0.135 

-0.493 

-0.185 

0.043 

0.060 

[M/H] 

-3.603 

0.109 

-0.433 

0.039 

0.027 

0.049 

[a/M] 

-4.360 

0.060 

-0.848 

-0.096 

0.013 

0.048 


Therefore, the ASPCAP abundances are not strictly 
differential to those of the Sun. Given the way the line 
lists were determined, the existence of the temperature 
trends in the raw abundances, and the lack of external 
calibration (apart from that for [M/H]), APOGEE de¬ 
fines its own relative abundance scale. We consider the 
accuracy of this scale by several external comparisons, 
including observations of the solar spectrum and of Arc- 
turus, observations of stars with independently-measured 
abundances, and consideration of the locus in paramters 
and abundances of the entire APOGEE sample. 

5.5.1. Solar and Arcturus abundances 

Observations of the solar spectrum, obtained by ob¬ 
serving the asteroid Vesta, and of Arcturus were made 
using the NMSU lm feed to the APOGEE spectrograph. 
These spectra are of particular interest because their 
abundances set the standard for many abundance stud¬ 
ies, they have high resolution spectral atlases, and be¬ 
cause the combination of the abundances and high reso¬ 
lution spectra were used to tune the APOGEE line list. 

However, the internal calibration that depends on the 
ASPCAP temperature was derived from the cluster sam¬ 
ple which is significantly cooler than the Sun. Addition¬ 
ally, the calibration was derived from giants, not dwarfs, 
so the applicability of the calibration to the solar abun¬ 
dances is debatable, and this complicates the solar com¬ 
parison. Table [7] presents the solar abundances derived 
from the Vesta spectrum, showing both the uncalibrated 
abundances and the calibrated abundances determined 
by applying the internal calibration correction appropri¬ 
ate for 5250 K, which is the high temperature end of the 
calibration sample. 

In general, the observed uncalibrated solar abundances 
are within the expected errors, except for magnesium, 
sulfur, and vanadium. Applying the calibration (again, 
outside of its derived range) slightly improves some of 
the abundances, but makes the magnesium, silicon, man¬ 
ganese, and vanadium abundances worse. Systematic off¬ 
sets from the solar abundance all appear to be within 0.1 
apart from sulfur and vanadium, and perhaps silicon (if 
one considers the calibrated value). 

Table [8] presents the ASPCAP Arcturus abundances, 
along with several different comparison values: the abun- 


Table 7 

Solar parameters and abundances from Vesta 


Element Uncalibrated ASPCAP Calibrated ASPCAP 


Ted 


5728.± 91.5 

5812.± 91.5 

log g 


4.54± 0.11 

4.49± 0.11 


M/H[ 

0.03± 0.03 

0.11± 0.03 


a/M 


0.00± 0.01 

0.02± 0.01 


C/H 


0.06± 0.07 

0.06± 0.07 


N/H 


0.10± 0.06 

0.10± 0.06 


O/H 


0.03± 0.01 

0.02± 0.01 


Na/H] 

0.03± 0.03 

0.03± 0.03 


Mg/H] 

0.07± 0.03 

0.11± 0.03 


Al/H] 

-0.03± 0.01 

0.00± 0.01 


Si/Hl 

0.03± 0.03 

0.12± 0.03 


S/HI 


0.26± 0.04 

0.22± 0.04 


K/H] 

-0.07± 0.03 

-0.04± 0.03 


Ca/H] 

0.02± 0.07 

-0.03± 0.07 


Ti/H] 

-0.07± 0.04 

0.05± 0.04 


V/H] 

0.11± 0.05 

0.24± 0.05 


Mn/H] 

0.04± 0.03 

0.19± 0.03 


Fe/H] 

0.01± 0.07 

0.05± 0.07 


Ni/H] 

0.01± 0.15 

0.10± 0.15 


dances that were adopted to tune the line list, a set of 
abundances derived from a high-resolution near-IR FTS 
spectrum by Smith et al. (2013) using a similar (but not 
identical) linelist to the current APOGEE linelist, and a 
set of abundances d erived fro m atomic transitions in the 
optical spectrum by Ramirez & Allende Prieto (2011). 

Assessing the quality ol the comparison is challenging 
because there is not a consensus in the literature for all 
Arcturus abundances. Note that we do not necessarily 
expect the ASPCAP results to match the abundances 
that were adopted to tune the line list, because the line 
list tuning involved making compromises between the so¬ 
lar spectrum and the Arcturu s spectrum (with the latter 
receiving lower weight, Shetro ne et al.|2015 ) and only al¬ 
lowed atomic transition probabilities to vary within lim¬ 
its. 

Several different potential issues can be noted. For 
carbon, the near-IR abundances agree, b ut differ signifi¬ 
cantly fr om the Ramire z fe Al lende Prieto (2011) value; 
however, Sneden et al. "(2014) find a value much closer 
to the near-IR values. For oxygen, the ASPCAP value 
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Table 8 

Arcturus parameters and abundances 


Element 

ASPCAP 

Adopted 

Smith 

Ramirez 

ASPCAP [X/Fe] 

Adopted [X/Fe] 

Smith [X/Fe] 

Ramirez [X/Fe] 

T e ff 


4295.93± 91.47 


4275.00 

4286.00 





log 9 


1.71± 0.11 


2.10 

1.66 





M/H] 


-0.64± 0.03 


-0.40 






a/ M] 


0.20± 0.02 


0.40 






C /H 


-0.43± 0.04 

-0.43 

-0.43 

-0.09 

0.15 

0.09 

0.04 

0.43 

N /H 


-0.52:1 0.08 

-0.14 

-0.14 


0.05 

0.38 

0.33 


O /H 


-0.32:1 0.03 

-0.02 

-0.02 

-0.02 

0.26 

0.50 

0.45 

0.50 

Na/H] 

-0.61± 0.12 

-0.52 


-0.41 

-0.04 

0.00 


0.11 

Mg/Hl 

-0.40± 0.03 

-0.12 

-0.38 

-0.15 

0.18 

0.40 

0.09 

0.37 

Al/H 


-0.32=t 0.07 

-0.22 

-0.21 

-0.18 

0.26 

0.30 

0.26 

0.34 

Si/Hl 


-0.33± 0.04 

-0.12 

-0.39 

-0.19 

0.25 

0.40 

0.08 

0.33 

S /HI 


-0.41± 0.04 

-0.12 



0.17 

0.40 



K /H] 

-0.54=t 0.06 

-0.52 

-0.29 

-0.32 

0.04 

0.00 

0.18 

0.20 

Ca/Hl 

-0.51± 0.04 

-0.42 

-0.47 

-0.41 

0.07 

0.10 

0.00 

0.11 

Ti/H 


-0.52:1 0.05 

-0.12 

-0.31 

-0.25 

0.06 

0.40 

0.16 

0.27 

V /H] 

-0.77± 0.13 

-0.52 

-0.39 

-0.32 

-0.20 

0.00 

0.08 

0.20 

Mn/Hl 

-0.59:1 0.04 

-0.52 

-0.53 

-0.73 

-0.02 

0.00 

-0.06 

-0.21 

Fe/H 


-0.58:1 0.04 

-0.52 

-0.47 

-0.52 

0.00 

0.00 

0.00 

0.00 

[Ni/H] 


-0.51± 0.04 

-0.52 

-0.46 

-0.46 

0.07 

0.00 

0.01 

0.06 


differs significantly from either the alternate near-IR or 
the optical abundance, and seems quite low given that 
Arcturus has generally been found to be enhanced in a- 
element abundances. This may partly be understood be¬ 
cause oxygen is primarily derived from OH li nes, which 
are quite sensitive to temperature (see Smith et al. [20131 
and the raw Arcuturus ASPCAP temperature (4207K) is 
lower than the literature values (here is a case where one 
might have gotten better results if the calibrated tem¬ 
peratures were used for the abundances). For other a 
elements (Mg, Si, Ca), t he ASPC AP values agree fairly 
well with the |Smith et al.| (|20 1 3[) va l ues, b ut are lower 
than the Ramirez & Aliende Prieto (2011) values, sug¬ 
gesting less overall a enhancement. 

5.5.2. Other individual star comparisons 

In addition to Arcturus and Vesta, we obtained NMSU 
lm + APOGEE observations of several other bright stars 
(<5 Oph, (3 And, /r Leo) for which Smith et al. (2013) have 
derived abundances from high-resolution near HR spectra 
obtained with an FTS spectrograph and using the AS¬ 
PCAP line list. We also have observations of a number 
of nearby st ars with publi s hed abundances from optica l 
ana lyses by ReddyetaL (2006), |Reddy et al. (2003), 


and|Hensby et al. (|2014|); because these stars are bright, 


observations were also all obtained with the NMSU lm 
+ APOGEE. However, these are mostly warmer stars 
with higher surface gravity and fall outside the range 
of most APOGEE main sample stars and thus outside 
the range where we have internally calibrated the abun¬ 
dances. We have also compiled a set of abunda nces of 


et al. 

2015 

tor more details), to which we can compare 

the APOGEE/ASPCAP abundances. 


Meszaros et al. 

(2015) have derived abundances from 


members. In addition, [Cunha et al. (|2015|) have derived 
Na, O and Fe abundances from manual analysis (i.e., do¬ 
ing synthesis around individual lines) of APOGEE spec¬ 
tra for stars in the cluster NGC 6791, so these provide an 
opportunity to test the APOGEE/ASPCAP abundances 
against those derived by more traditional methods from 


the same spectra. However, these independent analyses 
of the APOGEE spectra were done using a slightly differ¬ 
ent version of the APOGEE linelist, so exact agreement 
is not expected (but the changes should be small, apart 
from K). 

Results for all of these samples are shown in Figure [7J 
[8j[9| and[Io|for stellar parameters ( T e g and logg), CNO, 
a-element, and other elements, respectively. In these 
plots, uncalibrated parameters/abundances are shown in 
the left panels, and calibrated parameters/abundances 
are shown on the right, both plotted as a function of 
effective temperature. Since, as discussed above, stars 
with higher surface gravity are not calibrated, most of 
the warmer stars are missing for the calibrated parame¬ 
ters/abundances. In all plots, the points are color-coded 
by the derived overall metallicity, with blue representing 
metal-poor stars and red metal-rich ones. Different sym¬ 
bols represent different comparison samples: solid circles 
are the results from the FTS spectra, solid triangles from 
the optical analyses of solar neighborhood stars. Open 
triangles show the independent analysis for NGC 6791, 
and the small dots show the results for the independent 
analysis of globular cluster stars. Although it is challeng¬ 
ing to interpret all of the comparisons because the com¬ 
parison samples are heterogeneous and each have their 
own associated uncertainties, a few potential issues can 
be noticed. 

From Figure [TJ the ASPCAP carbon abundances for 
giants seem to be in reasonable agreement with litera¬ 
ture values. ASPCAP nitrogen abundances seem to be 
significantly lower, and oxygen abundances moderately 
lower. 

The a element abundances in Figure [8] suggest that Mg 
compares fairly well, that Si is a bit high in the ASPCAP 
abundances, and that Ca and Ti are a bit low. Tita¬ 
nium results may indicate a trend with temperature; as 
discussed below, there is other indication that ASPCAP 
titanium abundances may have significant issues that are 
not currently fully understood. 

For the metal-poor globular cluster stars, the ASPCAP 
temperatures, surface gravities, and abundances appear 
offset from the literature values and, in fact, even from 
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Figure 7. Comparison with literature spectroscopic effective tem¬ 
peratures and surface gravities. Left panel is uncalibrated AS- 
PCAP results; right is for calibrated results. Differences are 
ASPCAP minus r eference values. Large circles are abundances 
from FTS spectra l|Smith et al. 20 l".T|l. filled triangles from optical 
spectra and analys is, open tri angles fr om independent analysis of 
APOGEE spectra |[Cunha et al. 2015J. 

the values measured from the independent analyses of 
the APOGEE spectra. We suspect that this is related to 
non-standard abundance patterns in some of the globular 
cluster stars, e.g., oxygen abundances that do not scale 
with other a elements. The assumption that a elements 
vary in lockstep made during the initial stellar parameter 
determination may lead to stellar parameters that are 
systematically offset, which is evident in the comparison 
of effective temperatures for the globular cluster stars. 
The systematic offsets with the stellar parameters may 
then lead to systematic offsets in some of the individual 
abundances. In addition, since features such as CN and 
CO can be very weak in metal-poor stars, the adopted 
scheme of allowing [C/M] and [N/M] to vary with the 
stellar parameters may cause problems. 

Overall, these comparisons suggest that there may 
be non-negligble systematic offsets betweeen ASPCAP 
abundances and literature abundances for some ele¬ 
ments. Such offsets are perhaps not uncommon in abun¬ 
dance analysis where derived values have not been ad¬ 
justed relative to those obtained from a source of similar 
stellar parameters (e.g., the Sun for analysis of warmer 
dwarfs). Nonetheless, they do indicate that, in an exter¬ 
nal sense, the accuracy of ASPCAP abundances may be 
poorer (0.1-0.2 dex) than the internal precision. 

5.6. Sample characteristics 

In this section, we provide information about the global 
distribution of properties across the entire APOGEE 
sample, to help to indicate the quality of derived param¬ 
eters and abundances across the large parameter space 
of the APOGEE stars. 

5.6.1. Quality of the fits 

For each stellar spectrum that is fit by ASPCAP, 
FERRE returns a y 2 of the fit. The y 2 distribution for 
the entire calibrated APOGEE DR12 sample is plotted 
in Figure EH It is clear that the quality of the fits are 
poorer for cooler stars. This is perhaps not unexpected 
as the number and depth of spectral features generally 
increase at cooler temperatures, so any problems, e.g., 
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Figure 8. Same as Figure [7| for CNO abundances. 
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Figure 9. Same as Figure [q for a-element abundances (Mg, Si, 
S, Ca, Ti. 
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Figure 10. Same as Figure [7] for Na, Al, K, V, Mn, Fe, and Ni. 
with the line list, will be amplified at cooler tempera¬ 
tures. 

These results suggest that there is probably more un¬ 
certainty in the ASPCAP results at cooler temperatures. 
For each star, if the y 2 value is larger than the typical 
value at the temper ature o f the star, a bit is set in the 
ASPCAPFLAG (see f 7.1.31: the dashed and solid lines in 
Figure [TT] show the levels above which the CHI2_WARN 
and CHI2_BAD bits are set, respectively. 

5.6.2. Stellar parameters 

One way of testing the quality of the derived stellar 
parameters is to see to what extent we reco ver the ex¬ 
pected trends on an HR diagram. Figure [12] presents an 
HR diagram from both the uncalibrated APOGEE DR12 
sample (top) and the calibrated sample (bottom), with 
points color-coded by their derived metallicity, [M/H], 
We also overlay isochrones for a 4 Gyr old population 
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Figure 11. \ 2 distribution for full APOGEE giant sample. Lines 
show the limits above which the CHI2_WARN (dashed) and 
C'HI2_BAD (solid) bits are set. 


(Bressan et al.|[2012i). 

The uncalibrated results (top panel) show that the 
dwarfs, while clearly separated from the giants, have stel¬ 
lar parameters that are offset from the stellar models. In 
part, this is probably due to the exclusion of a stellar ro¬ 
tation dimension in the synthetic grid, but since not all 
stars are expected to have significant rotation, there are 
probably additional effects as well. While the ASPCAP 
parameters are of sufficient quality to separate dwarfs 
from giants, these issues, along with the lack of calibra¬ 
tors available for stars of high surface gravity, lead us 
to recommend that the parameters and abundances for 
dwarfs be used only with caution, and no attempt has 
been made to calibrate them. 

The bottom panel of Figure |T2] shows the calibrated 
giant sample. The calibrated parameters show the ex¬ 
pected change in the locus of points as a function of 
metallicity. This agreement, which is obtained without 
imposing any assumptions about the relation between 
effective temperature, surface gravity, and metallicity, is 
quite encouraging. We note that some scatter in the lo¬ 
cus at fixed metallicity is expected if there exists a range 
of ages at given metallicity; this full sample includes stars 
from across the Milky Way, so this is not expected to be a 
coeval sample. In particular, the stars with lower metal¬ 
licity may be older than the 4 Gyr old isochrones shown 
here; older isochrones would be in better agreement with 
the derived ASPCAP parameters. 

The problems with the surface gravities for red clump 
stars are visible, as stellar evolution models predict that 
red clump stars should have near constant surface grav¬ 
ity, while the ASPCAP sequence slopes to higher surface 
gravity at higher temperature. The ASPCAP gravities 
for RC stars are also higher than expected from the stel¬ 
lar models (as was also s een in the comparison to the 
asteroseismic gravities in (5.21. 
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Figure 12. Effective temperature - surf ace gr avity - me tallici ty distributions for uncalibrated full APOGEE sample (top) and for calibrated 
giant sample (bottom). Isochrones from Bressan et al.| ( |2012| ) for a 4 Gyr population are overlaid. 
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Figure 13. Metallicity - [a/M] - temperature relation for cali¬ 
brated APOGEE giant sample. 

Figure [l3| shows the [a/M] vs [M/H] relation for the 
calibrateogiant sample, color coded by temperature. 
The expected trends towards higher [a/M] at lower 
metallicity is apparent. However, the cooler stars (blue 
points) appear to deviate from the locus of stars of other 
temperatures. While it is possible that there are some 
astropliysical effects contributing here (cooler stars are 
more likely to be at larger distances), a more likely inter¬ 
pretation is that the [a/M] results (and probably other 
abundances) are less reliable at the cooler end of the AS- 
PCAP temperature range, i.e. T < 4000 K. This is also 
the temperature range where the quality of the fits gets 
significantly worse. As a result, users should exercise 
caution when using results for the cooler stars. 


5.6.3. Stellar abundances 


Figure 14 presents plots of [X/Fe] vs. [Fe/H] for all 15 
elements for a subsample of stars with high S/N (> 200); 
the scatter in these plots increases for lower S/N. These 
are presented to assess to what degree some of the ex¬ 
pected trends are observed; a more extensive discussion 
of these plots in terms of chemical evolution is beyond 
the scope of this paper. Points in the plots are color- 
coded by effective temperature; black points are those 
for which the empirical uncertainties ([5.4.2) are larger 
than 0.1 dex. 

Previous studies from near the solar neighborhood 
have shown, in general, enhanced [a/Fe] ratios at lower 
metallicities, while other elements tend to have less vari¬ 
ation. The APOGEE results for the a elements O, Mg, 
Si, S, and Ca show trends of increasing abundance rel¬ 
ative to Fe at lower metallicity, as expected, although 
the behavior of Mg at the lowest metallicities indicates 
some possible issues in that regime. However, Ti does 
not indicate any significant trend. Given that such a 
trend has been clearly o bserved for titan ium in the so¬ 
lar neighborhood (e.g., Bensby et al. 2014), this, along 


with the individual star comparison discussed previously, 
suggests that there is some issue that may be affecting 
the reliability of the ASPCAP Ti abundance. Possible 
issues with the ASPCAP derived Ti abundances will be 
investigated further in future work. 

While the shape of the a element locii are gener¬ 
ally as expected, the location of the suggest possible 
problems with the overall accuracy (zeropoint) of the 
APOGEE/ASPCAP abundances. The Si and Su abun¬ 
dances ratios ([Si/Fe] and [S/Fe]) are greater than zero at 
solar metallicity by 0.1-0.2 dex, while the Ca abundance 
([Ca/Fe]) is around -0.1. As noted above, some offsets 
are perhaps not unexpected since we do not reference our 
abundances to those of any well-calibrated object. 

Results for Na and V show large scatter; features from 
these elements are weak and most likely the scatter re¬ 
sults from limited precision of the APOGEE measure¬ 
ments. 

The points in Figure [14] are color-coded by the derived 
effective temperature, so that trends with temperature 
may be discerned. In general, it seems that temperature 
trends are relatively small, with the possible exception 
of those for Ti, Mn, and Ni. However, given that cooler 
giants are more luminous, stars of different temperatures 
likely sample different regions of the Milky Way, so ap¬ 
parent abundance trends with temperature do not nec¬ 
essarily represent an issue with the abundance determi¬ 
nations, as the chemical relations may vary across the 
Milky Way. 

Overall, these results support our previous assessment 
that the internal precision of APOGEE/ASPCAP abun¬ 
dances is fairly good, while the external accuracy is more 
uncertain. 


5.7. Persistence 

As noted in the description of the instrument, portions 
of two of the SDSS-III/APOGEE detectors have strong 
persistence that affects roughly one-third of the fibers 
over 40% of the APOGEE spectral range. Because the 
effect is complicated, not fully understood, and depends 
on the details of the exposure history, no attempt has 
yet been made to correct for it in SDSS-III/APOGEE 
analysis. 

We attempt to determine whether persistence has a 
significant effect on the derived stellar parameters and 
abundances by comparing results for objects whose spec¬ 
tra were recorded outside the persistence region in all 
visits against results for objects whose spectra were 
recorded within the highest persistence region in all vis¬ 
its; there are some intermediate cases where some visits 
are recorded in the persistence region, and some outside 
of it (since plates are generally replugged between visits), 
as well as cases where spectra fall on regions of interme¬ 
diate or lower persistence. 

Results are shown in Figure [l5j [l6j ITT] and [18] In all 
cases, we show only fainter stars with H > 11, which 
tend to be more affected by persistence than brighter 
stars. The left plots show the results for stars that do 
not fall in the persistence region for any of the visits, 
and the right plots show results for stars that fall in the 
high persistence region in all visits. Figure [15] suggests 
that persistence does not seem to impact the parame¬ 
ters dramatically. However, for the individual elemental 
abundances, persistence can have a noticeable effect. In 
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Figure 14. Abundance ratio plots for high S/N (>200) stars from full APOGEE giant sample. Points are color coded by effective 
temperature; black points are stars with large uncertainties (empirical uncertainties greater than 0.1 dex). 
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Figure 15. Persistence comparison for T e ^ vs loggfor stars with 
H > 11. Left panel has stars with all observations outside of 
persistence region, while right panel has stars with all observations 
inside of persistence region. Points are color-coded by temperature. 
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Figure 16. Persistence comparison for CNO abundances for stars 
with H > 11. Left panel has stars with all observations outside of 
persistence region, while right panel has stars with all observations 
inside of persistence region. Points are color-coded by temperature; 
black points have empirical uncertainties greater than 0.1 dex. 

particular, the relations for N, Mg, S, K, Ca, Ti, Mn, 
and Ni show significantly larger scatter or offsets for stars 
that fall in the high persistence region; given that some 
of these have features that are not within the persistence 
region, the stellar parameters must be affected at some 
level. 

As a result, users are cautioned about the use of some 
of the individual element abundances for stars that were 
observed within the persistence region. These stars can 
be identified using the star quality bitmask described be¬ 
low. 



[Fe/H] 


[Fe/H] 


Figure 17. As Figure [16| but for a-element abundances (Mg, Si, 
S, Ca, Ti). 

including low-level SQL access through the CasJobs in¬ 
terface at skyserver.sdss3.org/Casjobs. More infor¬ 
mation is available at skyserver.sdss3.org. 

The CAS provides access to several database tables: 

• apogeeStar — a table with basic information for 
each combined spectrum, including targeting flags, 
radial velocity information, and the STARFLAG 
and ANDFLAG data quality bitmasks for each in¬ 
dividual star; 

• aspcapStar — a table with the ASPCAP parame¬ 
ters, abundances, uncertainties, and bitmasks flag¬ 
ging possible conditions with the parameters and 
abundances; 

• apogeeVisit — a table with the radial velocity in¬ 
formation for each individual visit of each star; 


6. DATA PRODUCTS AND ACCESS 

APOGEE data can be accessed in several different 
ways. 

6.1. Catalog Archive Server (CAS) 

Catalog data are also stored in a database called the 
Catalog Archive Server (CAS), which can be accessed 
at http://skyserver.sdss3.org/drl2. The CAS con¬ 
tains a number of ways of interacting with the database, 


• apogeeObject — a table with all of the ancillary 
catalog/targeting information for the parent sam¬ 
ple of stars in each field from which the observed 
targets were drawn; 

• apogeePlate — a table with information about each 
plate observed; 

• apogeeDesign — a table that contains information 
about each plate design; 
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Figure 18. As Figure [To] but for Na, Al. K, V. Mn, and Ni. 

• apogeeField - TTtable with information about each 

field that was observed, which includes the LOCA- 
TION_ID numeric value that is used to index dif¬ 
ferent fields, and also the planned number of visits 
for each field. 

There are also some tables that link these tables, e.g., 
apogeeStar to apogeeVisit. 

This database only contains catalog information de¬ 
rived from the spectra, but not the spectra themselves. 

6.2. Science Archive Server (SAS) 

The data files created at multiple stages in the data re¬ 
duction and processing pipeline are available on the Sci¬ 
ence Archive Server (SAS), which is accessible through 
the web at http://data.sdss3.org/sas/drl2 

6.2.1. FITS tables 

Summary FITS tables are available that present com¬ 
pilations of data for all stars in the data release. These 
FITS tables are the original data source for the CAS 
database tables: 


• allStar — summary information on each individ¬ 
ual star observed in SDSS-III/APOGEE, which in¬ 
cludes mean barycentric radial velocity, ASPCAP 
parameters and abundances as derived from the 
summed spectra, and a compilation of ancillary tar¬ 
geting data (for DR12, allStar-v603.fits); 

• allVisit — summary information on each individual 
visit of each star in SDSS-III/APOGEE, which in¬ 
cludes the barycentric radial velocity of each visit, 
along with a compilation of ancillary targeting data 
(for DR12, allVisit-v603.fits). 

• allPlates — summary information for all observed 
plates, include the location and plate design infor¬ 
mation (for DR12, allPlates-v603.fits). 


Both allStar and allVisit files are FITS files with sev¬ 
eral header/data units (HDUs). The first HDU con¬ 
sists of the main data, with an entry for each star (all¬ 
Star HDU1) or visit (allVisit HDU1). The objects are 
sorted by increasing right ascension. The second HDU 
in both files contains a 360 element integer array giving 
the HDU1 array index of the first object larger than each 
degree of right ascension, which may be useful for a rapid 
search through the table for an object of known RA. The 
third HDU contains information about the order of array 
values in the PARAM (PARAM.SYMBOL) and ELEM 
(ELEMJ3YMBOL and ELEMLVALUE) arrays in HDU1; 
these are constant for all objects and thus do not need 
to be repeated for each object. 

The allStar file provides pointers into the allVisit file to 
identify the visits that went into the combined spectra. 

There are also several summary targeting FITS tables: 

• apogeeDesign gives targeting information for each 
plate design 


• apogeePlate gives targeting information for each 
plate (designed, location_id, hour angle, tempera¬ 
ture, epoch) 

• apogeeField gives targeting information for each 
field (name, location_id, coordinates, and expected 
number of visits) 


• apogeeObject gives all of the ancillary cata¬ 
log/targeting information for all stars in each field 
from which the observed targets were drawn 


Targeting terms are defined in Zasowski et al. (2013), 
particularly in the glossary of that paper. 


6.2.2. Spectral data 

The spectral data files created at multiple stages in the 
data reduction and processing pipeline are also available 
on SAS. Some of the data that might be of the most 
interest include: 


• apVisit files — FITS image files with individual 
visit spectra of each star on each MJD in which it 
was observed; 

• apStar files — FITS image files with combined 
(across multiple visits) spectra of each star ob¬ 
served with APOGEE, along with uncertainty and 
mask arrays, as well as the individual visit spectra 
resampled to the apStar rest wavelength scale; 
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• aspcapStar files — FITS image files for each star 
that give the continuum normalized spectra from 
which parameters and abundances are derived, the 
adopted uncertainties in the spectra, and the best¬ 
matching synthetic spectra, along with derived pa¬ 
rameters and abundances in the image header; 

• apField files — FITS tables for each field that in¬ 
clude information and the radial velocities for all 
stars in the field; 

• aspcapField files — FITS tables for each field 
that include the derived ASP CAP parameters and 
abundances for all stars in the field, along with the 
spectra and best-matching synthetic spectra. 

These files are described in more detail in {js| Note that 
the information in the apField and aspcapField files is 
compiled including all fields in the allStar file discussed 
above. 

In addition to these main products, the SAS contains 
all of the data from APOGEE, from the raw data cubes 
to the final processed spectra and derived parameters, 
including all of the intermediate data products. 

Table [9] provides a rough guide to the overall data 
structure in the SAS. All of these data products are 
documented through the SDSS datamodel, available at 
http://data.sdss3.org/datamodel 

6.3. Web interface to the SAS 

A web application provides an interface to the SAS 
through a front-end database called the SAS database 
(SASDB). This provides a tool by which individual spec¬ 
tra can be interactively inspected without having to 
download the data, and also provides an interface by 
which you can search for data to download either for in¬ 
dividual objects, or for a bulk search of multiple objects. 
This can be found at http: //data. sdss3. org 

7. USING APOGEE DATA 

7.1. Bitmasks, data quality, and target classes 

SDSS uses bitmasks to flag stars with quality issues 
that may affect the accuracy of the derived parameters 
and abundances, as well as to provide information on tar¬ 
geting. A bitmask is an integer in which individual bits 
(corresponding to powers of two in a decimal representa¬ 
tion) can be set to indicate a variety of conditions. Users 
are advised to routinely check the contents of these. 


7.1.1. Targeting bitmasks 

There are different target classes for objects that 
appear in APOGEE data releases. While the bulk of the 
targets were chosen for the “main” survey (i.e., stars se¬ 
lected with a simple color selection), there are a number 
of other targe t classes, as described in detail in |ZasowsE| 
et al. (20131. These include special targets such as 
the cluster targets used for the calibration, targets for 
ancillary science programs, and observations obtained 
with the NMSU lm feed to the APOGEE spectrograph. 
The targeting class for each target is documented in a 
pair of targeting bitmasks, APOGEE_TARGETl and 
APOGEE_ TARGET2, whe r e the bit definitions are 
defined in Zasowski et al. (|2013| and documented at 
http://www.sdss.org/dr12/algorithms/bitmasks/ 


#AP0GEE_TARGET1 and http://www.sdss.org/drl2/ 
algorithms/bitmasks#AP0GEE_TARGET2. We provide 
character “translations” of the targeting bitmasks in the 
APOGEE.TARGFLAGS column in summary data files. 

Note that the main survey targets are not defined 
by a single bit in the targeting bitmasks; instead, 
they can be identified as targets that have one of 
the APOGEEJ3HORT, APOGEEJNTERMEDIATE, or 
APOGEEJLONG bits set. For convenience, we have in¬ 
cluded a separate bitmask, EXTRATARG, in the sum¬ 
mary files and in the CAS tables. This bitmask has a 
value of zero for all main survey targets, but has dif¬ 
ferent bits set if the target is either a star selected for 
science outside of the main survey (bit 0), commission¬ 
ing data (bit 1), a telluric star (bit 2), or a star observed 
with the NMSU lm (bit 3), see http://www.sdss.org/ 
dr12/algorithms/bitmasks#APOGEE_EXTRATARG. 

There are some objects that may appear more than 
once in the summary files or the CAS tables. Many of 
these are stars that were observed during the commis¬ 
sioning period. These observations are kept separate be¬ 
cause the quality of the the data is poorer (worse LSF). 
Since these data are not homogeneous with the main sur¬ 
vey data, they are not analyzed with ASPCAP; they are 
included because they still provide useful radial veloc¬ 
ity information and may reveal spectral variations over 
time (e.g., double-line spectroscopic binaries or emission 
line stars). Apart from this, there are ~ 1600 stars that 
appear more than once with survey data. This occurs be¬ 
cause these objects were targeted in different, but over¬ 
lapping fields. In some cases, this was intentional (in 
particular, for some of the ancillary science programs), 
but for others the targets were inadvertently selected in¬ 
dependently in two different overlapping fields. For all 
duplicate targets, we choose a “primary” observation as 
the observation with the highest S/N, and all other ob¬ 
servations of this target have a bit (bit 4) set in the EX¬ 
TRATARG bitmask. 

The most homogeneous sample are those stars with 
EXTRATARG=0. This will give only main survey tar¬ 
gets, excluding all special targets, commissioning data, 
lm objects, telluric stars, and duplicate targets. 

7.1.2. STARFLAG bitmask 

Each visit spectrum has an associated STARFLAG 
bitmask, w hich flags a number of conditions as given 
in Table 1 10| and at http://www.sdss.org/algorithms/ 
drl2/bitmasks/#AP0GEE_STARFLAG. If BAD.PIXELS, 
VERY_BRIGHT_NEIGHBOR, or LOW_SNR is set for a 
visit spectrum, that spectrum is deemed bad and not 
used to obtain the final combined spectrum. The PER- 
SISTTIIGH, PERSISTJV1ED, and PERSISTJLOW bits 
flag stars that may be significantly affected by the per¬ 
sistence on the blue detector; they are set if 20% or more 
pixels are in the persisence region, with the three dif¬ 
ferent levels representing different amplitude of the per¬ 
sistence affect (note that PERSIST JLOW is still signifi¬ 
cant persistence: stars outside the persistence region al¬ 
together do not have any of the these bits set). 

The spectra that are combined from the individual visit 
spectra have STARFLAG set to the logical OR of the 
STARFLAG values of the individual visits. In addition, 
there is another ANDFLAG that is set to the logical 
AND of the flags from the individual visits. 
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Table 9 

SAS directory structure 


Directory 


contents 


APRED.VERS 

APRED_VERS/red/MJD 

APRED_VERS / apo25m/PL ATE/M JD 
APRED_VERS/apolm/PLATE/MJD 
APRED.VERS / APSTAR.VERS 

APRED.VERS / APSTAR.VERS / apo25m/LOC ATION.ID 
APRED.VERS/APSTAR_VERS/apolm/PROGRAM_ID 
APRED.VERS /APSTAR.VERS / ASPC AP.VERS 

APRED.VERS /APSTAR.VERS /ASPC AP.VERS/location_id 
APRED.VERS /APSTAR.VERS /ASPC AP.VERS/RESULTS.VERS 

APRED.VERS /APSTAR.VERS / ASPC AP.VERS/RESULTS.VERS/location_id 


top level directory for a given visit reduction version 
reduced data for individual exposures, 
ap2D and aplD 

reduced visit data for 2.5m data: apVisit files 
reduced visit data for lm data: apVisit files 
top level for a given star combination version 
combined, resampled spectra: apStar files for 2.5m data 
combined, resampled spectra: apStar files for lm data 
top level directory for a given ASPCAP 
version / configuration 
FERRE files for a given location 

top level directory for a given RESULTS (calibration) 
version: allStar and all Visit summary tables 
ASPCAP results in FITS files for given field: 
aspcapStar and aspcapField files 




Table 10 

STARFLAG bitmask 


Condition 

Bit 


Description 


BAD-PIXELS 

0 

COMMISSIONING 

1 

BRIGHT-NEIGHBOR 

2 

VERY-BRIGHT-NEIGHBOR 

3 

LOW-SNR 

4 


5-8 

PERSIST-HIGH 

9 

PERSIST-MED 

10 

PERSIST-LOW 

11 

PERSIST-JUMP-POS 

12 

PERSIST-JUMP-NEG 

13 


14-15 

SUSPECT-RV-COMBINATION 

16 

SUSPECT-B ROAD JUNES 

17 


18-31 


Spectrum has many bad pixels (>40%): BAD 

Commissioning data (MJD<55761), non-standard configuration, poor LSF: WARN 
Spectrum has neighbor on detector more than 10 times brighter: WARN 
Spectrum has neighbor on detector more than 100 times brighter: BAD 
Spectrum has low S/N (S/N<5): BAD 
Currently unused 

Spectrum has significant number (>20%) of pixels in high persistence region: WARN 
Spectrum has significant number (>20%) of pixels in medium persistence region: WARN 
Spectrum has significant number (>20%) of pixels in low persistence region: WARN 
Spectrum shows obvious positive jump in blue chip: WARN 
Spectrum shows obvious negative jump in blue chip: WARN 
Currently unused 

WARNING: RVs from synthetic template differ significantly 
from those from combined template 

WARNING: cross-correlation peak with template significantly 
broader than autocorrelation of template 
Currently unused 


7.1.3. ASPCAP bitmasks 


Several bitmasks are used in the ASPCAP processing. 
There is an overall ASPCAPFLAG, which flags various 
conditions as listed in Table 11 In particular, bits are 
set if any of the parameters are near the spectral library 
grid edges, if the ASPCAP temperature differs signifi¬ 
cantly from that expected from the observed colors, if 
the spectrum shows signs of significant stellar rotation, 
or if the spectrum does not meet the survey requirement 
for S/N. 

Of particular note is the STAR-BAD bit, which is set if 
any of other bits (TEFF-BAD, LOGG-BAD, CHI2_BAD, 
COLORTE_BAD, ROTATION-BAD, SN_BAD) to de¬ 
note that the ASPCAP results are very likely not to be 
reliable. Stars with this bit set should not be used for 
parameter/abundance analysis. Note that some bits in 
ASPCAPFLAG are informational/warning only, so re¬ 
stricting a sample to ASPCAPFLAG=0 is not an appro¬ 
priate filter to get good data. 

In addition, there are PARAMFLAG and ELEMFLAG 
arrays for each of the stellar parameters and individ¬ 
ual element abundances that provide additional informa¬ 
tion about why a given parameter may be denoted BAD 


in ASPCAPFLAG or that flag potential issues with the 
abu ndance determination for individual elements. Table 
[l2| gives the definition of the PARAMFLAG bits. 

7.2. Radial velocities 

Several measurements of radial velocities are provided, 
as discussed in the following subsections. 


7.2.1. Systemic radial velocities 

The average systemic velocity for each APOGEE star 
is given in VHELIO-AVG; estimated uncertainties are 
given in VERR (weighted error), and VERR-MED (me¬ 
dian error of individual visit RVs). The scatter obtained 
from multiple visits (with number of visits given in NVIS- 
ITS) is given in VSCATTER. Note that objects with 
VSCATTER > 0.5 km/s are potentially RV-variable ob¬ 
jects, base d on the typical RV accuracy (see Nidever| 


et al. 


2015). These velocities were determined by deriv- 


mg relative visit RVs by iterative cross-correlation of each 
visit against the combined spectrum for the object, with 
an absolute scale established after convergence by cross¬ 
correlation of combined spectra with the best-matching 
synthetic spectrum in a RV template grid. This has the 
advantage of having no template mismatch for the rela- 
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Table 11 

ASPCAPFLAG bitmask 


Condition Bit 


Description 


TEFF.WARN 

0 

LOGG.WARN 

1 

VMICRO .WARN 

2 

M.H.WARN 

3 

ALPHA.M.WARN 

4 

C.M.WARN 

5 

N.M.WARN 

6 

STAR.WARN 

7 

CHI2.WARN 

8 

COLORTE.WARN 

9 

ROTATION.WARN 

10 

SN.WARN 

11 

12-15 

TEFF.BAD 

16 

LOGG.BAD 

17 

VMICRO .BAD 

18 

M.H.BAD 

19 

ALPHA _M .BAD 

20 

C.M.BAD 

21 

N.M.BAD 

22 

STAR.BAD 

23 

CHI2.BAD 

24 

COLORTE.BAD 

25 

ROTATION.BAD 

26 

SN.BAD 

27 

28-30 

NO.ASPCAP .RESULT 

31 


WARNING on effective temperature (see PARAMFLAG[0] for details) 

WARNING on log g (see PARAMFLAG[1] for details) 

WARNING on vmicro (see PARAMFLAG[2] for details) 

WARNING on [M/H] (see PARAMFLAG [3] for details) 

WARNING on [a/M] (see PARAMFL AG 4] for details) 

WARNING on [C/M] (see PARAMFL AG 5] for details) 

WARNING on [N/M] (see PARAMFL AG [6] for details) 

WARNING overall for star: set if any of TEFF, 

LOGG, CHI2, COLORTE, ROTATION, SN warn are set 
high x 2 (> 2*median at ASPCAP temperature) (WARN) 
effective temperature more than 1000K from photometric 
temperature for dereddened color (WARN) 

Spectrum has broad lines, with possible bad effects: ratio of FWHM of cross-correlation 
of spectrum with best RV template to FWHM of auto-correlation of template >1.5 (WARN) 
S/N < 70 (WARN) 

Not currently used 

BAD effective temperature (see PARAMFL AG [0] for details) 

BAD log g (see PARAMFLAG [1] for details) 

BAD vmicro (see PARAMFLAG[2] for details) 

BAD ' ' 

BAD 
BAD 
BAD 


3] for details) 

4] for details) 

5] for details) 

6] for details) 

BAD overall for star: set if any of TEFF, LOGG, CHI2, COLORTE, ROTATION, 

SN error are set, or any parameter is near grid edge (GRIDEDGE_BAD is set in any PARAMFLAG) 


M/H] (see PARAMFLAG 
a/M] (see PARAMFLAG' 
C/M] (see PARAMFLAG 
N/M] (see PARAMFLAG 


high x 2 (> 5*median at ASPCAP temperature) (BAD) 

effective temperature more than 2000K from photometric temperature for dereddened color (BAD) 
Spectrum has broad lines, with possible bad effects: ratio of FWHM of cross-correlation 
of spectrum with best RV template to FWHM to auto-correlation of template > 2 (BAD) 

S/N < 50 (BAD) 

Not currently used 

True for all commissioning data, and for spectra that are marked as bad 


Table 12 

PARAMFLAG bitmask 


Condition 

Bit 

Description 

GRIDEDGE.BAD 

0 

Parameter within 1/8 grid spacing of grid edge 

CALRAN GE.BAD 

1 

Parameter outside valid range of calibration determination 

OTHER.BAD 

2 

Other error condition 


3-7 

Not currently used 

GRIDEDGE.WARN 

8 

Parameter within 1/2 grid spacing of grid edge 

CALRAN GE.WARN 

9 

Parameter in possibly unreliable range of calibration determination 

OTHER.WARN 

10 

Other warning condition 


11-15 

Not currently used 

PARAM.FIXED 

16 

Parameter set at fixed value, not fit 


17-31 

Not currently used 


tive RVs (although a disadvantage is that the combined 
template is not noiseless). The parameters of the best¬ 
matching grid point are given in RV.TEFF, RVJLOGG, 
and RV_FEH, but these are not to be confused with the 
much better determined ASPCAP parameters; no inter¬ 
polation is performed in the RV grid. 

An alternative systemic velocity is given in SYN- 
THVHELIO.AVG, with uncertainties in SYNTHVERR 
and SYNTH VERR.MED and scatter in SYNTH VSCAT- 
TER. These are determined by cross-correlation of each 
visit with the RV template grid, and then subsequent 
combination. 

The scatter between the two different systemic RV 
measurements is given in SYNTHSCATTER. Larger val¬ 
ues here suggest possible issues with the derived RVs. For 
such objects, the SUSPECT_RV.COMBINATION bit in 


the STARFLAG bitmask is set. 

7.2.2. Individual visit radial velocities 

The individual visit radial velocities are stored in the 
visit tables (allVisit FITS table in the SAS, and the al- 
IVisit table in the CAS). The visit RV information de¬ 
rived from cross-correlating each visit against the com¬ 
bined spectra are stored with labels VREL (the measured 
RV, uncorrected for Earth motion), VRELERR (uncer¬ 
tainty in VERR), and VHELIO (measured RV corrected 
for barycentric motion; we note the slight misnomer in 
the VHELIO label since these are strictly barycentric, 
not heliocentric velocities). The visit RVs derived from 
cross-correlating each visit against the best-matching 
template are stored in SYNTH VREL, SYNTH VREL¬ 
ERR, and SYNTH VHELIO. 
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An initial estimate of the radial velocity for each for 
each visit is made by cross-correlation of the individual 
visit spectra with the same coarse grid of synthetic tem¬ 
plates mentioned above. These are stored in the visit 
tables as ESTVREL, ESTVRELERR, and ESTVHELIO. 


7.3. Stellar parameters 

As described previously, each APOGEE spectrum is 
matched to a spectral grid covering several dimensions 
- T e ff, log g, v micro , [M/H], [C/M], [N/M], [a/M]— to 
derive a set of stellar parameters that are subsequently 
used for the individual element abundance fits. 

Calibrations to T e e , logg, [M/H], and [a/M], as de¬ 
scribed above, are applied to the raw FERRE out¬ 
put; for [M/H], both internal (temperature) and ex¬ 
ternal calibration is applied. These quantities are 
stored in the PARAM array in the allStar file (HDU1); 
the order of the parameters in the array is given 
in the PARAM_SYMBOL array in HDU3 of the all¬ 
Star file. In the CAS, the parameters are given 
in individual columns PARAMJTEFF, PARAMJLOGG, 
PARAMJLOG VMICRO, PARAMJVLH, PARAMALM, 
PARAMJNLM, and PARAM_ALPHA_M in the aspcap- 
Star table. 

The calibrated effective temperatures and surface grav¬ 
ities are duplicated in columns TEFF and LOGG in 
both the allStar file and CAS tables. These are the rec¬ 
ommended values for these quantities. The allStar file 
also duplicates the parameter-level [M/H] and [a/M] in 
PARAMJVLH and PARAM.ALPHAJM. 

The raw, uncalibrated FERRE output is stored in 
an FPARAM array in the allStar file, and in the 
CAS in columns FPARAMJTEFF, FPARAM JLOGG, 
FPARAM JM_H, FPARAM_CJM, FPARAMJNJM, and 
FPARAM_ALPHA_M in the aspcapStar table. 

There is a distinction between the parameter-level 
[M/H] and [a/M] values and the individual elemental 
abundances discussed below. The parameter-level quan¬ 
tities are fits to the entire spectrum, and might weight 
different elements differently over a range of temper¬ 
atures. However, they might also carry higher accu¬ 
racy as averages over multiple elements and more wave¬ 
lengths. Furthermore, the parameter-level [M/H] is the 
only abundance that is externally calibrated, based on 
the literature cluster metallicities discussed previously. 

Stars with unreliable parameters are flagged by having 
the STARJ3AD bit in the ASPCAPFLAG bitmask set. 
If this bit is set, then other bits in ASPCAPFLAG can 
be used to determine the reason for the measurements 
being flagged as unreliable. The STARJWARN bit can 
be consulted to determine if derived parameters may be 
suspect based on less stringent criteria. 

As discussed previously, while the ASPCAP results 
can easily discriminate dwarfs from giants, the results 
for dwarfs are significantly less secure and lack an ex¬ 
tensive calibration sample. As a result, no calibratrated 
parameters are provided for dwarfs. 

Based on the quality of the fits and external compar¬ 
isons, we also note that ASPCAP results at cooler tem¬ 
peratures, T e s< 4000 K may be more uncertain. 


7.4. Stellar abundances 

For the 15 indi vidual elemental abundances, the in¬ 
ternal calibration (5.4.1) of abundance with temperature 


is applied, and the results are stored in the ELEM ar¬ 
ray in the allStar file (HDU1); the order of the elements 
in the array is by atomic number and is given in the 
ELEM.SYMBOL array in HDU3 of the allStar file. Since 
carbon, nitrogen, and the a elements are fit by varying 
the [C/M], [N/M], and [a/M] dimensions in the spec¬ 
tral library, respectively, the abundances from FERRE 
for these elements are relative to the overall scaled-solar 
metal content, [M/H], while the abundances for other 
elements are relative to hydrogen since they are fit by 
varying [M/H]. These differences in the array contents 
are documented in the ELEMJVALUE array of HDU3, 
which explicitly gives what values are presented in the 


arrays. 

For convenience, we convert all calibrated abundances 
into [X/H] values (by adding [M/H] to the carbon, ni¬ 
trogen and a element abundances) and store these in 
named quantities: C_H, N_H, 0_H, Na_H, Mg_H, Al_H, 
Si_H, S_H, K_H, Ca_H, Ti_H, V_H, Mn_H, Fe_H, Ni_H in 
the allStar file and also as individual colum ns in t he as¬ 
pcapStar table in the CAS. As discussed in [ 5.4.3 these 
abundances do not have any external calibration applied 
to them. 

The uncalibrated abundances are stored in an 
FELEM array in the allStar file, and in individual 
columns FELEM_C_M, FELEM_N_M, FELEM_0_M, 
FELEM_Na_H, FELEMJMgJM, FELEM_AlJH, 
FELEMJSiJM, FELEMJ3JM, FELEMJKJH, 

FELEM_Ca_M, FELEM_Ti_M, FELEM_V_H, 
FELEMJMnJH, FELEM_Fe_H, FELEM_Ni_H in the 
CAS table. Note that these are not all converted into 
[X/H]; as raw values, they preserve the way in which 
they were derived. 


The ASPCAPFLAG should be consulted to confirm 
that there are no significant issues with the stellar pa¬ 
rameters. In addition, before using the abundance of 
any particular element, it is important to confirm that 
the derived value is not near an abundance grid edge. 
This is achieved by consulting the appropriate ELEM- 
FLAG bitmask, looking to see if the GRIDEDGEJBAD 
bit (value within 1/8 of a grid spacing of edge) or the 
GRIDEDGEAVARN bit (value within 1/2 of a grid spac¬ 
ing of edge) is set. This is especially critical when looking 
at abundance ratios, where the value of the ratio may not 
immediately indicate a problem with one or both of the 
abundances. 

Similarly to the stellar parameters, abundances for 
dwarfs are more uncertain, and calibrated abundances 
for dwarfs are not provided. 

The calibrated giant abundances have a eff ectiv e tem¬ 
perature correction applied, as described in §5.4| Since 
this metallicity-independent liner correction is perhaps 
unlikely to represent the full behavior of the abundances, 
a sample that covers a wide temperature range may have 
larger possibility for temperature-dependent systemat- 
ics than one that is restricted to a smaller temperature 
range. In addition, the calibration sample only extends 
down to T e ff~ 3800 K, so the behavior at cooler temper¬ 
atures is less well known; the current results apply the 
offset at 3800 K for all cooler stars, but set the CAL- 
RANGEJWARN bit in the ELEMFLAG bitmask. 


Results for some elements may be less reliable than 
for others. In particular, the titanium abundances show 
large scatter compared with previous measurements for 
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a calibration subsample, and do not show the expected 
behavior with overall metallicity. Other elements like 
sodium and vanadium have relatively few and weak fea¬ 
tures, and thus have large r unc ertainties. 

In addition, as noted in §5.7[ individual element abun¬ 
dances for objects whose spectra fell in the the persis¬ 
tence region of the blue detector may be more uncer¬ 
tain. These objects can be identified by those having 
one of the PERSIST .HIGH, PERSISTJV1ED, or PER- 
SISTJLOW bits set in the ANDFLAG bitmask (for ob¬ 
jects falling in persistence region in all visits), or those 
with one of the persistence bits set in the STARFLAG 
bitmask (for objects falling in persistence region in any 
of the visits). 

The empirical error estimates derived from simple fits 
to the observed scatter of abundances in clusters as a 
functi on of temperature, metallicity and signal-to-noise 
(see 15.4.1) are stored in the ELEM_ERR array (and 
correspondingly named columns X_H_ERR in the CAS). 
Note that the scatter used for these estimates is the scat¬ 
ter in the raw FERRE abundance measurements, i.e., in 
[X/M] for C, N and the a elements, but in [X/H] for Na, 
Al, K, V, Mn, Fe, and Ni. The FELEM_ERR array (and 
corresponding column names in the CAS) gives the raw 
FERRE errors. 

For uncertainties in abundance ratios, it is not neces¬ 
sarily appropriate to add the uncertainties of the individ¬ 
ual abundances in quadrature, as the scatter may not all 
result from random errors. If there are systematic errors 
contributing to the scatter, it is possible that different 
elements will be offset in the same way, such that abun¬ 
dance ratios could be more robust. Further analysis on 
this issue is planned for the future. 


7.5. Selection effects 

As with all surveys, there are selection effects in the 
APOGEE catalog. Several effects to consider are: 


• main sample selection: only a fraction of the stars 
in each field are selected for observation, with a 
wel l-defined distribution i n magnitude and color — 
see 


Zasowski et al. (2013) for details. 


• ancillary and calibration targets: in addition to the 
main survey targets, there are calibration and an¬ 
cillary program targets that are not selected by 
the usual algorithm. While the non-main survey 
targets can be recognized from the bitmasks (e.g., 
EXTRATARG=0), there is a subtle effect that any 
targets that are part of an ancillary program are 
rejected from being selected as a part of the main 
sample, so it is possible that some object types are 
actually underrepresented in the main sample. 


• abundance biases at warmer temperatures: the 
sample selection generally uses a (dereddened) blue 
color cut of (J — K )o > 0.5. While this was de¬ 
signed to include metal-poor luminous giants, it 
does bias against lower metallicity stars on the 
lower giant branch. Furthermore, the abundance 
measurements may become less certain at warmer 
temperatures because features are weaker. 


• abundance biases at cooler temperatures: the cur¬ 
rent spectral grids are limited to T e // > 3500 K 


and logg> 0, so no cooler and lower surface grav¬ 
ity stars appear with good ASPCAP parameters. 
This biases against the most metal-rich stars, in 
particular, for the most luminous giants as well as 
bright AGB stars. Furthermore, the abundance 
measurements may become less certain at cooler 
(T e // < 4000 K) temperatures because of the com¬ 
plexity of the spectra and current limitations in our 
ability to model it well. 

7.6. Caveats 

There are issues that are discovered in public data 
releases after the point when things can be modified. 
The online SDSS/APOGEE documentation attempts to 
track these at http://www.sdss3.org/drl2/irspec/ 
caveats. At the time of the data release, known caveats 
include: 


• Uncertainties given as [X/H] are actually computed 
from [X/M] scatter for C, N and a elements. 

• Character translations for target flags (in 

APOGEE.TARGFLAGS) are missing in the 
ap Visit and allVisit files for four ancillary 
programs (APOGEE_RV_MONITOR_IC348, 

APOGEE_RV_MONITOR_KEPLER, 
APOGEE_GES_CALIBRATE, 
APOGEE_BULGE_RV.VERIFY); they are correct 
in the apStar and allStar file 

• Some targets observed in multiple fields and appear 
as duplicates (as flagged with the EXTRATARG 
flag discussed previously). 

• Some targets may have been selected indepen¬ 
dently for different programs within different visits 
to the same held. As a result, there are stars con¬ 
structed from multiple visits for which a target hag 
bit may be set in the combined spectra, but not in 
all of the visit spectra that were used to construct 
it. 


• The APOGEE.ANCILLARY bit is not 
set for APOGEEJCEPLERJ3EISMO and 
APOGEE.KEPLERJHOST targets. 

• For a few objects, slightly incorrect object/star 
names were used during the reduction. 


8. USING APOGEE SPECTRA 
8.1. Raw data 

Raw data hies are saved as apR-[abc]-ID8.apz frames, 
where [abc] refers to each of the three chips, and the 
ID8 is an eight digit number that uniquely identihes each 
APOGEE exposure. The hrst four digits give the number 
of days since 2010 December 31, i.e. MJD-55562. The 
last four digits are a running number of the frames taken 
during that day. 

Data from the three detectors are stored in separate 
FITS hies, with the longest wavelengths at the lowest 
columns in chip a and the shortest wavelengths at the 
highest columns in chi p c. Each hie is com pressed using 
the FPACK algorithm (Seaman et al. 2010); the resulting 
images are given as an initial image plus a set of differ¬ 
ences between adjacent up-the-ramp reads. We highlight 
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Table 13 

Contents of ap Visit files 


Table 14 

Contents of apStar files 


HDU contents 


0 main header with target information, starflag 

1 flux 

2 uncertainty 

3 pixel bitmask 

4 wavelength 

5 sky 

6 sky uncertainty 

7 telluric correction spectrum 

8 telluric correction uncertainty 

9 wavelength coefficients 

10 LSF coefficients 


HDU contents 


0 main header with target, RV, version information 

1 resampled spectra 

2 uncertainties in spectra 

3 pixel mask for spectra 

4 sky spectrum 

5 sky uncertainty 

6 telluric spectrum 

7 telluric uncertainty 

8 LSF coefficients 

9 RV information, including template 
cross-correlation and autocorrelation 


the fact that some reconstruction is necessary to give the 
raw data cubes by giving the apR files a .apz extension 
name (rather tha n .fits), even though they are valid FITS 
files. See|Nidever et al. (20151 for more details. 


8.2. Visit spectra 

Individual visit spectra of each star on each MJD 
in which it was observed are saved in ap Visit- VERS- 
PLATE-MJD-FIBER.fits files. VERS represents the 
data release version of the reduction routines, which is 
r5 for the DR12 release. PLATE is the SDSS plate iden¬ 
tification that corresponds to a physical set of holes into 
which fibers are plugged. MJD represents the modified 
Julian date of the observations, and FIBER represents 
which fiber was used to obtain the spectrum. The FIBER 
number is defined such that the spectrum at the high¬ 
est row numbers has FIBER=1, while th ose at the lovw 
est ro w numbers have FIBER=300 (see Majewski et al. 
20151. The association of the fiber number with an object 


name is provided by a “plugmap” file; this association is 
used to populate header information with an object name 
and coordinates. 

The apVisit spectra combine the multiple dithered 
exposures taken on a given night. They are on the 
native APOGEE pixel scale, but with half-pixel sam¬ 
pling made possible by the dithered exposures. Mul¬ 
tiple extensions give uncertainty in the spectra, a bit- 
mask flag for individual pixel status (http: //www. sdss. 
org/dr12/algorithms/bitmasks/#APOGEE_PIXMASK ), a 
wavelength array, the sky and telluric spectra used to 
correct the object spectrum, and some other ancillary 
information. The content of the files is summarized in 
Table [H 

There will be some bad pixels in any given spectrum 
that result from bad pixels on the APOGEE detectors. 
Such pixels are flagged in the pixel bitmask that is stored 
in HDU3 of the apVisit files, which encodes the reason 
why any given pixel is flagged as unreliable. The bitmask 
definitions for these are given at http: //www. sdss. org/ 
drl2/algorithms/bitmasks/#APOGEE_PIXMASK. The 
bitmask also encodes locations where data may be more 
suspect, e.g., pixels in the strong persistence regions of 
the detect ors, and/or pi xels i n the location of the Littrow 
ghost (see [Wilso n et al.||2015 |. 

While correction is made for telluric absorption and for 
sky emission, these corrections (especially for sky emis¬ 
sion) are generally imperfect. As a result, users can ex¬ 
pect to see some residual sky emission in the visit spectra. 


Locations of strong sky emission and telluric absorption 
are flagged in the pixel bitmask. 

8.3. Combined spectra 

FITS image files with combined (across multiple visits) 
spectra of each star observed with APOGEE are saved in 
apStar-VERS-APOGEE_ID.fits files. Here, VERS repre¬ 
sents the reduction version (again, r5 for the DR12 re¬ 
lease) and APOGEEJD represents an object identifier, 
which is usually the 2MASS object name with a ’2M’ 
prefix; for the few objects that are not in the 2MASS 
catalogs, a similar style APOGEEJD was created from 
the coordinates, and a ’AP’ prefix was prepended. 

The apStar files represent the main survey output 
spectra. These FITS files have multiple HDUs: they 
include the combined spectra and individual resampled 
visit spectra, uncertainty arrays corresponding to these, 
a pixel mask array that flags pixels with various con¬ 
ditions, a wavelengt h ar ray, and other information, as 
summarized in Table IT4l 

The combined spectrum is included along with the in¬ 
dividual visit spectra shifted to rest wavelength in HDU1. 
This is a 2D image with wavelength running along the 
column axis, and with NVISITS+2 rows that have two 
versions of the combined spectra in the first two rows, 
followed by NVISITS rows of resampled individual-visit 
spectra. The two combinations are generally very sim¬ 
ilar; they differ i n the w eighting schem e used to do the 


combination (see Nidever et al. 

2015 

). The ASPCAP 

analysis ol these spectra uses the 

first version ot tfie com- 


bined spectra. Bad pixels in any of the individual visit 
spectra are discarded during the combination. 

All spectra are put on a common wavelength scale with 
a fixed dispersion in log A: apStar spectra all have 8575 
pixels with log Xi = 4.179+6. x 10 -6 i, where i is the pixel 
number (with the first pixel having * = 0). This corre¬ 
sponds to a pixel sampling of 4.145 km/s and provides 
roughly 3 pixels per resolution element. 

Because of the resampling, the uncertainties are cor¬ 
related between adjacent pixels; the covariance informa¬ 
tion is not currently being tracked. If data with non- 
correlated errors are required, users are advised to use 
the apVisit files, but these have non-constant dispersion 
and, even here, there will be small covariances in the 
dither-combined images because the dithers are not per¬ 
fect. 

Because spectra are obtained with visits taken over a 
period of time, telluric absorption and sky emission fea- 
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Table 15 

Contents of aspcapStar files 


HDU contents 


0 Normalized APOGEE spectrum 

1 Uncertainties in normalized spectrum 

2 Best library match as derived by ASPCAP/FERRE 

3 FITS table with parameters of best fit spectrum 


Table 16 

Contents of aspcapField files 


HDU contents 


0 FITS table with parameters and abundances for all 

objects in field 

1 FITS table with normalized spectra, uncertainty array, 

and best-fit spectra 

2 FITS table giving the order of the parameters and 

abundances in the arrays of HDUO 


tures generally do not fall at the same rest wavelength 
in the stellar spectra because of the changing barycen- 
tric correction. This means that systematic errors in the 
correction of these features are generally reduced com¬ 
pared to the individual visit spectra, but are spread over 
a larger region of the stellar spectra. The pixel bitmask 
in HDU3 of the apStar files track pixels that were af¬ 
fected in the individual visit spectra. 


8.4. ASPCAP spectra and best fits 

The aspcapStar files are FITS image files that give 
the continuum normalized spectra from which param¬ 
eters and abundances are derived. These are given on 
the apStar wavelength scale, but only include the spe¬ 
cific wavelength ranges that are matched to the syn¬ 
thetic grids. Multiple extensions give the best-matching 
synthetic spectra f rom which the stellar parameters are 
adopted (see Table 151. 


8.5. ASPCAP field results 

The aspcapField files are FITS tables for each field 
that include the derived ASPCAP parameters and abun¬ 
dances for all stars in the field. A second extension in¬ 
cludes a FITS table with all of the normalized spectra 
and the best-matching synthetic spectra for all stars in 
the field. Because these files include spectra for hundreds 
of stars, they are moderately large (see Table 16). 

9. CONCLUSION 


We have described the data from the SDSS- 
III/APOGE E surv ey as released in the Data Release 12 
(|Alam et al.||2015|). 

Stellar parameters and abundances have been cali¬ 
brated and compared with the literature using obser¬ 
vations of stellar clusters and of stars with previously 
measured abundances. The internal precision of the 
APOGEE abundances is typcially 0.05-0.1 dex, judging 
from the internal scatter of abundances within clusters. 
External accuracy of the abundances is challenging to 
assess, but may be good only to 0.1-0.2 dex. Stellar 
parameters and abundances for cooler stars (T e ff<4000) 
may be more uncertain. 


The APOGEE data set is rich and we expect that 
it will enable significant developments in Milky Way 
science. We have described the methods by which 
APOGEE data can be accessed, and attempt to present 
information that will provide information about its reli¬ 
ability and limitations. 

The automatic determination of stellar parameters and 
abundances presents a challenging problem, and we do 
not claim to be currently doing the most optimal job. 
Work is ongoing with a number of issues, including: 

• implementation of a correction for persistence in 
the SDSS-III/APOGEE data; 

• improvement of sky emission correction; 

• improvement of radial velocities, especially for hot¬ 
ter stars and spectra with lower S/N; 

• continued revision of line lists; 

• incorporation of rotation where needed; 

• extension to cooler temps usin g the 

_ MARCS/T urbospectrum code (see Zamora 

| et ahl|2015t 

• investigation of individual element grids, rather 
than using the element family grids to solve for 
individual element abundances; 

• investigation of correlated errors; 

• improvments in handling of non-detections/upper 
limits; 

• investigation of issues with derivation of surface 
gravity. 

Nevertheless, the DR12 results have proven to be useful 
and of a quality to conduct a wide variety of science. 

As the study of Milky Way stars continues with SDSS- 
IV/APOGEE-2, there will be future data releases, and 
we anticipate that at least some of these will include 
improved re-analysis of SDSS-III/APOGEE data. 
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